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ON THE GORTLER INSTABILITY IN HYPERSONIC FLOWS: 
SUTHERLAND LAW FLUIDS AND REAL GAS EFFECTS 1 

Yibin B. Fu, Philip Hall, and Nicholas D. Blackaby 
Department of Mathematics 
Oxford Road 
University of Manchester 
Manchester M13 9PL 
U.K. 

Abstract 

The Gortler vortex instability mechanism in a hypersonic boundary layer on a curved wall 
is investigated in this paper. Our aim is to clarify the precise roles of the effects of boundary 
layer growth, wall cooling and gas dissociation in the determination of stability properties. 
We first assume that the fluid is an ideal gas with viscosity given by Sutherland’s law. It 
is shown that when the free stream Mach number M is large, the boundary layer divides 
into two sublayers: a wall layer of 0(M 3 / 2 ) thickness over which the basic state temperature 
is 0(M 2 ) and a temperature adjustment layer of 0(1) thickness over which the basic state 
temperature decreases monotonically to its free stream value. Gortler vortices which have 
wavelength comparable with the boundary layer thickness (i.e. have local wavenumber of 
order M ~ 3 / 2 ) are referred to as wall modes. We show that their downstream evolution is 
governed by a set of parabolic partial differential equations and that they have the usual 
features of Gortler vortices in incompressible boundary layers. As the local wavenumber 
increases, the neutral Gortler number decreases and the centre of vortex activity moves 
towards the temperature adjustment layer. Gortler vortices with wavenumber of order one 
or larger must necessarily be trapped in the temperature adjustment layer and it is this 
mode which is the most dangerous. For this mode, we find that the leading order term in 
the Gortler number expansion is independent of the wavenumber and is due to the curvature 
of the basic state. This term is also the asymptotic limit of the neutral Gortler numbers 
of the wall mode. To determine the higher order correction terms in the Gortler number 
expansion, we have to distinguish between two wall curvature cases. When the wall curvature 
is proportional to (2x) -3 / 2 where x is the streamwise variable, the Mach number M can be 
scaled out of the problem and we show that in the 0(1) wavenumber regime, Gortler vortices 
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Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665.Additional support 
was provided by USAF under Grant AFOSR89-0042 and SERC. 
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are again governed by a set of parabolic partial differential equations and therefore the higher 
order correction terms in the Gortler number expansion are not uniquely determined and 
are strongly dependent on nonparallel effects. In the large wavenumber limit, however, 
nonparallel effects become of second order; Gortler vortices evolve downstream in a quasi- 
parallel manner and the Gortler number expansion has its first three terms independent of 
nonparallel effects. In the more general case when the wall curvature is not proportional to 
(2a:) -3 / 2 , the effect of the curvature of the basic state persists in the downstream development 
of Gortler vortices; non-parallel effects are important over a larger range of wavenumbers and 
they become of second order only when the wavenumber is of order higher than 0(M X / 4 ). 
In the latter case the Gortler number expansion has the first two terms independent of 
nonparallel effects; the first term being due to the curvature of the basic state and the 
second term due to viscous effects. The second term becomes comparable with the first term 
when the wavenumber reaches the order M 3 / 8 , in which case another correction term can 
also be found independently of nonparallel effects. Next we investigate real gas effects by 
assuming that the fluid is an ideal dissociating gas. We find that both gas dissociation and 
wall cooling are destabilizing for the mode trapped in the temperature adjustment layer, 
but for the wall mode trapped near the wall the effect of gas dissociation can be either 
destabilizing or stabilizing. 
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1 Introduction 


This paper is an extension of our previous paper Hall and Fu (1989) on the linear development 
of Gortler vortices and the reader is referred to that paper for a more detailed review of the 
relevant literature. In that paper, we assumed that the fluid was an ideal gas with the viscosity 
given by Chapman’s law. It was found that when a hypersonic boundary layer first loses 
stability to Gortler vortices, the vortices are necessarily trapped in the logarithmically thin 
temperature adjustment layer over which the temperature of the basic flow changes rapidly 
to its free stream value. In other words, the mode trapped in the temperature adjustment 
layer has a smaller Gortler number than any other mode. As a consequence of this region of 
vortex activity being thin (which leads us to consider Gortler vortices of small wavelength), 
the perturbation equations governing the downstream development of the vortices reduce to 
ordinary differential equations within the order of approximation considered if the appropriate 
“fast” streamwise dependence of the instability is built into the disturbance flow structure. 
Thus the non-uniqueness of the neutral stability curve associated with incompressible Gortler 
vortices disappears at high Mach numbers and a unique neutral curve with distinct left and 
right branches is obtained. 

However, a real fluid has its viscosity given by the more complicated Sutherland’s law. 
Although in most of the previous investigations on compressible boundary layers Chapman’s 
viscosity law has been adopted as an approximation to Sutherland’s law, such an approximation 
is poor for hypersonic flows in which the fluid temperature varies significantly across the 
boundary layer. Thus it is of interest to investigate how our previous results are modified 
if the more realistic Sutherland’s law is adopted. This is one of the problems which we are 
addressing in the present paper. 

The other problems which we consider are the effects of gas dissociation and wall cooling 
on the flow stability. For a hypersonic boundary layer, the temperature near the wall is 
typically of order 0(M 2 ) where M is the free stream Mach number, and gas dissociation must 
necessarily take place. Also, in practical situations, walls can not possibly withstand such high 
temperatures and they must be cooled. Thus it is also of special interest to clarify the precise 
roles of these two mechanisms in the stability properties of hypersonic boundary layers. 

The major difference between Gortler vortices in incompressible and hypersonic flows is that 
the presence of the temperature adjustment layer at the edge of a hypersonic boundary layer 
where the basic temperature field decreases rapidly to its free-stream value enables hypersonic 
Gortler vortices to be concentrated well away from the wall. In the incompressible case we 
know from the work of Hall (1982,1983), and Denier, Hall and Seddougui (1990) that at order 
one Gortler numbers unstable Gortler vortices are not localized within the basic boundary 
layer. At higher Gortler number the most dangerous Gortler vortices have a wavelength small 
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compared to the boundary layer thickness and are trapped near the wall. At order one Mach 
numbers this situation does not change significantly and the nonparallel problem has been 
discussed by Wadey (1990) and Spall and Malik (1989). In the latter two investigations the 
nonparallel equations were solved numerically following the approach of Hall (1983), the main 
result obtained was that the neutral position or growth rate of a Gbrtler vortex is a function 
of it’s upstream history. However the numerical calculations of Wadey (1990) suggest that as 
the Mach number increases the position where an unstable Gortler vortex locates itself moves 
towards the edge of the boundary layer. That result is consistent with what we shall find in 
this paper. 

The present paper is limited only to the linear regime of vortex growth ; nonlinear aspects of 
incompressible or low Mach number Gbrtler vortex growth have been discussed by Hall (1988), 
and Wadey (1990). For a detailed account of the nonlinear regime the reader is referred to 
the review paper by Hall (1990). Before going on to discuss the work presented here we 
also note that hypersonic boundary layers are also susceptible to instabilities not induced by 
streamline curvature. Thus, for example, Cowley and Hall (1990), Blackaby, Cowley and Hall 
(1990), Smith and Brown (1990) have discussed the role of Rayleigh or Tollmien-Schlichting 
wave instabilities in hypersonic boundary layers. Thus any nonlinear investigation of Gbrtler 
vortices at hypersonic speeds must allow for the possible interaction of the vortices with other 
finite amplitude instability mechanisms. 

This paper is organized as follows. In §2 and in the first part of §3 we discuss the basic 
state. Here we discuss the significant changes in the basic state which occur when Sutherland’s 
law is used instead of the Chapman law. In particular, the logarithmically thin temperature 
adjustment layer found for Chapman fluids is now replaced by a more complex adjustment layer 
of 0(1) thickness. We shall see later that this difference will strongly affect the downstream 
evolution properties of Gbrtler vortices. 

In §3 we formulate the linear stability problem for a hypersonic boundary layer and the 
linear perturbation equations are obtained in the usual manner by superimposing a Gbrtler 
vortex structure on the basic state and linearizing the Navier-Stokes equations. When Chap- 
man’s law is used, these perturbation equations can be reduced to a set of ordinary differential 
equations after the wavelength of the vortices is scaled by the thickness of the logarithmically 
thin adjustment layer. However, here we show that because the temperature adjustment layer 
is of 0(1) thickness when Sutherland’s law is used, non-parallel effects are more pronounced 
and their effects are different for different wall curvatures. To be more specific, when the 
wall curvature is proportional to (2z) -3 / 2 , the curvature of the basic state only gives rise to 
an 0(M 3 / 2 ) wavenumber independent term in the Gbrtler number expansion and its effect is 
not present in the downstream development of Gbrtler vortices in the neighbourhood of the 
neutral position; whilst in the more general case when the wall curvature is not proportional 
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to (2a:) -3 / 2 , the curvature of the basic state not only gives rise to an 0(M 3 / 2 ) wavenumber in- 
dependent term in the Gortler number expansion, but it also affects the downstream evolution 
of Gortler vortices in the neighbourhood of the neutral position and thus affects the determi- 
nation of other higher order correction terms to the Gortler number expansion. Sections 4 and 
5 are respectively devoted to the discussion of these two cases. 

In section 4, we show that in the special curvature case k(x ) = (2x) -3 / 2 , the Mach number 
can be scaled out of the problem; in the 0(1) wavenumber regime, the perturbation equa- 
tions are partial differential equations and they have to be solved numerically by a marching 
procedure. We present our numerical results which show that neutral curves depend crucially 
on what initial perturbations we impose and where we impose them. In the large wavenum- 
ber limit, nonparallel effects are negligible and a simple asymptotic expression is obtained for 
the Gortler number in terms of the wavenumber. In section 5, we show that because of the 
persistent effect of the curvature of the basic state, non-parallel effects are important over a 
larger range of wavenumbers and they become negligible only when the wavenumber is of order 
larger than 0(M 1 / 4 ). Thus for wavenumbers of order 0(M a ) with a < 1/4, the perturbation 
equations which govern the downstream evolution of Gortler vortices are partial differential 
equations and the situation is similar to the 0(1) wavenumber case discussed in section 4. 
When a > 1/4, nonparallel effects are not so pronounced and the Gortler number expansion 
has the first two terms independent of nonparallel effects; the first term due to the curvature 
of the basic state and the second term due to viscous effects. The second term becomes com- 
parable with the first term when the wavenumber reaches the order 0(M 3 / 8 ), in which case 
another correction term can also be found independently of nonparallel effects. 

To complete our stability analysis, we devote section 6 to the wall mode which has wave- 
length comparable with the boundary layer thickness. This mode is nonparallel and neutral 
curves have to be obtained by solving a set of partial differential equations. Our numerical 
results show that neutral curves, although non-unique, all decrease monotonically with the 
wavenumber and tend to a constant value in the large wavenumber limit, thus matching in the 
large wavenumber limit with the mode trapped in the temperature adjustment layer. 

In section 7, we investigate real gas effects and wall cooling effects. We assume that the 
fluid is an ideal dissociating gas. After dissociation has taken place, the fluid becomes a gas 
mixture. We first determine the constitutive properties of the gets mixture and then show how 
our previous results for ideal gases are modified when gas dissociation is taken into account. 
We show for the mode trapped in the temperature adjustment layer that the leading order 
Gortler number is decreased by both gas dissociation and wall cooling and thus we conclude 
that both these mechanisms are destabilizing. For the wall mode, neutral curves are not 
unique and so we cannot draw any general conclusion. For the case we consider, the neutral 
curves corresponding to the two models intersect, so the effect of gas dissociation can be either 
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destabilizing or destabilizing. Finally in tbe last section we give some further discussion. 


2 Basic state 

Consider a hypersonic boundary layer over a rigid wall of variable curvature (1 /A)k(x*/L), 
where L is a typical streamwise length scale and A is a lengthscale characterizing the radius of 
curvature of the wall. We choose a curvilinear coordinate system (x*,y*, z") with x* measuring 
distance along the wall, y* perpendicular to the wall and z* in the spanwise direction. The 
corresponding velocity components are denoted by (u*,v*,w*) and density, temperature and 
viscosity by p*,T* and p* respectively. The free stream values of these quantities will be 
signified by a subscript oo. We define a curvature parameter 6 by 


and consider the limit S — * 0 with the Reynolds number R defined by 

u t> L P*oo 


R = 


taken to be large so that the Gortler number 


P*oo 


( 2 . 1 ) 


( 2 . 2 ) 


G = 2 RV 2 6 


(2.3) 


is 0(1). In the following analysis, coordinates (x*,y*,z*) are scaled on (£*,iZ -1 / 2 £, R~ l l 2 L), 
the velocity (u*,v*,w*) is scaled on (v,l 0 ,R~ 1 ^ 2 ul 0 ,R~ 1 ^ 2 ul 0 ) and other quantities such as 
p*,T*, and p* are scaled on their free stream values with the only exception that the pressure p * 
is scaled on p^u*^ and the coefficient of heat conduction k* is scaled on p^. All dimensionless 
quantities will be denoted by the same letters without a superscript *. Then the Navier-Stokes 
equations are given by 


dp d . . 

* = °' 

Du dp d , du. d.du 
p ~Dt ~ ~~di + dy^dy^ + dz^dz^* 


Dv 1 , dp d 

>TH+2 a ™ 3 = - R ‘ir, + Ty 


' (\ 2 \ dV P 


+ 


d dv p 
dxp ^ dy ^ 


d . dv. d . dv. 
+ dy^dy^ + dz^dz^’ 


Dw dp 9 2 . 

, W = - Se 9l + ft{ (A -3' ,) 


dw s 


dvp 1 d dv p d dw d 
dx p J + dx fi ^ dz^ + dy^dy^ + dz^dz^’ 


!*,% = M7 - DM’K^) 2 + (£)>] + (1 - 1)M’[1 - 


dy‘ 


,du , 

' dz ‘ 


,dh s 

■dp' 


Dp 

Dt 


(2.4) 

(2.5) 


( 2 . 6 ) 

(2.7) 
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, 1 d ,,dT. 1 d n dT. , . 

+ a dy( k dy^ + a dz^ k dz^’ ^ 

lM 2 p = (1 + a)pT. (2.9) 

Here we have used a mixed notation in which (vi,V2,V3) is identified with (u,v,w) and 
( xi,X2,X3 ) with ( x,y,z ). Repeated suffices signify summation from 1 to 3. The functions 
A, k, Cp and h denote in turn the bulk viscosity, the coefficient of heat conduction, the specific 
heat at constant pressure and the enthalpy per unit mass. The constants 7 ,M and a are in 
turn the ratio of specific heats, the Mach number and the Prandtl number defined by 


7 = 




M = 


uz 


U, 


OO 
2 > 


Poo c poo 

° ~ ~k ’ 

A. 00 


c Poo 73 tTZg a‘ x 

where SR is a gas constant and Oqo = \Jl is the sound speed in the free stream. Finally, the 
function a in the equation of state (2.9) denotes the percentage by mass of the mixture which 
has been dissociated. Later in §7 we shall give the expression for a for a specific dissociation 
model used in our discussion. In equations (2. 5) -(2. 8), the operator D/Dt is the material 
derivative and it has the usual expression appropriate to a rectangular coordinate system. 
The basic state is given by 

( u , v, w) = ( u(x,y ), v(x,y), 0), T = T(x,y), 


p=p(x,y), p = fi(x,y). 


( 2 . 10 ) 


By substituting (2.10) into the governing equations (2.4)-(2.9) it is straightforward to obtain 
the reduced equations satisfied by the basic state. The reader is referred to the book by 
Stewartson (1964) for a detailed discussion of these basic state equations. If we define the 
Howarth-Dorodnitsyn variable y and a similarity variable r) by 

V- f pdy and 7/ = -|L, (2.11) 

Jo \/2x 

then the continuity equation is satisfied if u and v are written as 

« = /'(*?). v = + f'(r)) r -zdr}\. (2.12) 

V2a: p Jo P 

Here the functions f(r]) and T(t ]) must satisfy 


a" + m n y = o, 

-(pkT'y + CpfT' + £(7 - 1 )M 2 p(f") 2 = 0, 

l T 


(2.13) 

(2.14) 


if the i-momentum and energy equations are to be satisfied. These equations must then be 
solved such that /, f vanish at the wall, /', T = 1 at infinity and either T ' = 0 or T specified 
at the wall. The y-momentum equation gives 

!?=» 

ay 
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to leading order so that p = p(x). In our following analysis, we assume that there is no pressure 
gradient along the streamwise direction and therefore we can take p = constant. Equation (2.9) 
then reduces to 

[1 + a(T)]pT = 1. (2.15) 

Note that in obtaining (2.13) and (2.14) we have not made any constitutive assumptions so 
that they are valid for dissociated gases (to be discussed in §7) as well as for undissociated 
ideal gases (to be discussed in §3— §6). 

3 The perturbation equations 

We first assume that the fluid is an ideal (one component) gas undergoing no dissociation so 
that a = 0. Then we can assume that (i). the specific heats are constants; (ii). the coefficient 
of heat conduction is linearly related to the shear viscosity and (iii). the enthalpy h is given 
by h = CpT. These assumptions lead to the results 

k = p., Cp = 1, p=i. (3.1) 

(Note that all of these quantities have been non-dimensionalized). Then the basic equations 
(2.13) and (2.14) simplify to 

//" + (f/")' = 0, (3.2) 

+ ff + (7 - 1 )M 2 |(/") 2 = 0. (3.3) 

These two equations can then be solved if we make an constitutive assumption about the 
viscosity /2. In the previous paper, Hall and Fu (1989), we used Chapman’s viscosity law. Here 
we use Sutherland’s viscosity law, the dimensionless form of which is given by 

7 * 3/2 

^ = + ( 3 - 4 ) 

1+771 

where 771 is a constant. Equation (3.4) is exact in the sense that it is derivable from the kinetic 
theories of gases (see Chapman and Cowling (1970) for a discussion of its validity, also compare 
(3.4) with (7.9)). At high Mach numbers we know from the work of for example Freeman and 
Lam (1959) that the basic state splits up into two distinct regions. Near the wall a boundary 
layer forms in which the downstream velocity approaches it’s free-stream value of unity whilst 
the temperature decreases from it’s value at the wall and ultimately decays algebraically at 
the edge of the layer. In the next region this algebraic decay is taken up and the temperature 
then approaches exponentially the free-stream value of unity. 

As mentioned above an explicit analytical solution for the equations (3.2) and (3.3) is not 
possible. However, an asymptotic analysis in the large Mach number limit shows that the 
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boundary layer can be divided into two regions: an inner region 77 = and an outer 

region 77 = 0(1). In the inner region, we define Y,T and F(Y) by 

Y = M^ 2 rj, T = M~ 2 T, F(Y) = M 1 / 2 /. (3.5) 

Then equations (3.2) and (3.3) give 

p 11 

(l + m)(-=)' + FF" = 0, (3.6) 

Vt 

i ^ ( ^ ) ' + J ? T' + h-l )tl + » ) ^=° (3.7) 

to leading order. These equations are to be solved numerically subject to the conditions 

F(Q) = f'(0) = 0, f (00) = 0, F'{ 00) = 1, 

T\ 0) = 0 if the wall is thermally insulated, (3.8) 

T’(O) = nT w if the wall is under cooling, 

where T w is the wall temperature scaled on M 2 T oa when the wall is thermally insulated and n 
is the wall cooling coefficient. 

In Fig.l we have shown the results of our numerical integration of the wall layer equations 
(3.6) and (3.7). The temperature profiles are plotted for three values of the wall cooling 
coefficients: n=0.2, 0.6 and 1.0 and were calculated with 7 = 1.4, a — .72, m = .509. The 
asymptotic profile for large Y given by equation (3.10) is also plotted there for comparison. 

For large Y, equations (3.6) and (3.7) have the asymptotic solutions 


F - Y ^+(r_ /) )3/. +-■ 

(3-9) 

r 9(i + *) 2 1 

a 2 (Y - )S ) 4 ’ 

(3.10) 


where both P and D are to be determined by a numerical calculation. The numerical values 
of P corresponding to four values of the wall cooling coefficients are listed in Table 1 together 
with the values of .F"(0), T'(0) and T(0). 

Table 1 



n = 0.2 

n = 0.4 

n = 0.6 

n = 0.8 

n = 1 

f "( 0) 

0.1517 

0.1997 

0.2317 

0.2560 

0.2758 

T'( 0) 

1.8192 x 10- 2 

1.7976 x 10- 2 

1.3909 X 10- 2 

7.6872 x 10~ 2 

0 

n o) 

3.286 x 10 -2 

6.572 x 10- 2 

9.858 X 10- 2 

1.3144 x 10- 1 

1.643 X 10 

p 

3.1808 

2.8366 

2.6301 

2.4840 

2.3721 
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(3.11) 


The asymptotic expressions (3.9) and (3.10) imply that in the region 77 = 0(1), 


1 77 + 

T = T{ri) + • • • . 

On substituting (3.11) and (3.12) into (3.2) and (3.3), we obtain to leading order 


(1 + m ) 


iT^n +,r=o ' 


These two equations are to be solved numerically subject to the matching conditions 

„ s, . D - 9(1 + mf 1 

as 77 — ► 0 , /(*?)- ^ 77 , ^ ^ ^ + 


and the conditions at infinity 


/'( 00) = 0 , T(oo)=l. 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


The result from such a numerical calculation is shown in Fig.2. In this figure the asymptotic 
result is the one given by (3.15). 

It can be shown from (3.13) and (3.14) that 

")(^r <•■»> 

so that after (3.14) has been solved numerically, the function can be computed easily 

from this equation. Also, we note that whilst the solution of (3.14) is independent of the inner 
region solution and thus of the conditions at the wall, the function / is dependent on the inner 
region solutions through the matching constant D. 

We now assume, as in Hall and Fu (1989), that the flow is perturbed to spanwise periodic 
stationary vortex structure with constant wavenumber a. The linearized stability equations 
for these Gortler vortices are then found by linearizing the Navier-Stokes equations (2.4)-(2.9) 
about the basic state and retaining the leading order terms in the high Reynolds number limit. 
We obtain 


i{uU x + vUy) + (p.a 2 + b)U - ( fiUy) y + luyV 
- + ™y) + (A«v)v} T - JxUyTy = 0, 


(3.18) 
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l(v x + kuG)U + \pyU x - l -pU xy - p x ll y + l(uV x + vVy) + (pa 2 + b)V - |( pV y ) y 

1 12 4 

+ P V - [^(™* + VV y + \kGU 2 ) + -pU xy - -pyVv + ~(pVy)y + p X Uy]T 

2 4 2 1 

- pUy T x - [--pu x + -pv y )Ty + -pyiaW - -ia pWy = 0, (3.19) 

1 _ _ 1 2 

pxiaU + -piaU x + pyiaV + -piaV y - iaP - -p(u x + v y )iaT 

o o O 

- j(uw x + vWy) - lfia 2 W + (PWy)y = 0, (3.20) 

+ vf y )T - ^( u x + Vy)T - y^(uT x + VTy) 

1 1 _ w 

+ f(U x + Vy) - ^(T X U + TyV) + ia( T ) = 0, (3.21) 

if x u - 2(7 - l)M 2 pUyU y + IfyV + i(uT x + vTy) + ± a 2 T 

- {fS^* + vTy) + (7 - 1 )M 2 PU 2 V + ^(pTy)y^T - ^PTyTy - ^(pTy)y = 0. (3.22) 

Here p = dp/dT, whilst (U,V,W), P and T denote the vortex velocity field, pressure and 
temperature, respectively. Equations (3.18)— (3.22) differ from equations (2.11a-e) given in 
Hall and Fu (1989) only in that the bulk viscosity is taken to be zero here; that assumption is 
actually implied in that paper. 

It was shown by Hall (1982) that in the incompressible case the neutral curve for small 
wavelength vortices has G ~ a 4 and that the vortices are confined to a layer of depth a -1 / 2 
where the flow is locally most unstable. Hall and Malik (1989) extended this approach to the 
above system for M = 0(1) and wrote 

G = goa 4 + gia 3 + • • • . 


They found that the leading order growth rate a 2 6* has 6* given by 

.u 2 % 


** = £ + ( _ 
a K 2T 3 


v ««»%_ 

- ^ 2)90- 


(3.23) 


In the neutral case, S* = 0 and (3.23) then determines the neutral Gortler number go as a 
function of rj. The most unstable location rf is where go has its minimum. In Hall and Fu 
(1989), it is found that whei. Chapman’s law is used, rf* moves away from the wall as the Mach 
number increases. It is also found that the basic state temperature is 0(M 2 ) over most of 
the boundary layer and decreases rapidly to its free-stream value over a logarithmically thin 
adjustment layer sitting at the edge of the boundary layer. It is in this thin layer that rf 
lies and hence where go has it’s smallest order of magnitude in the large Mach number limit. 
Thus it is concluded that the thin temperature adjustment layer is most susceptible to Gortler 
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vortices. From the preceding discussion in this section we see that when Sutherland’s law is 
used, the logarithmically thin temperature adjustment layer corresponding to Chapman’s law 
is now replaced by a more complex temperature adjustment layer of 0(1) thickness. If we 
still use (3.23) with 6* = 0 to calculate the orders of go in the inner layer 77 = 0(M -1 / 2 ) and 
the outer temperature adjustment layer 77 = 0(1), we find that go = 0(M 15 / 2 ) in the former 
and go = 0(1) in the latter. Hence once again the temperature adjustment layer is most 
susceptible to Gortler vortices with wavenumber of order one or larger. It should be noted, 
however, the above conclusion is based upon a large wavenumber argument. In section 6 we 
shall show that the wall layer 77 = 0(Af -1 / 2 ) is actually of order M 3 / 2 thickness in terms of 
the physical variable y. Thus Gortler vortices with wavelength comparable with the boundary 
layer thickness must be trapped in the wall layer and have a = 0(M~ 3 / 2 ). It will be shown 
in section 6 that this wall mode has neutral Gortler number decreasing monotonically and 
has the centre of Gortler vortex activity moving towards the temperature adjustment layer as 
the wavenumber increases. Therefore, the minimum Gortler number corresponds to the mode 
trapped in the temperature adjustment layer and the latter is indeed the most dangerous mode 
when the whole range of wavenumbers are considered. It should also be noted that the result 
go = 0(1) for the temperature adjustment layer is obtained by taking the large Mach number 
limit of the 0(1) Mach number results. By doing so we have actually missed a term related 
to the curvature of the basic state which is not important for the case M = 0(1) and a 1, 
but is important in the large Mach number limit. As we shall show later on, the curvature 
of the basic state produces an effective Gortler number of order M 3 / 2 in the absence of wall 
curvature so that instability can not occur for G = 0(1). 

3.1 The strongly unstable inviscid mode 

Let us first confine our attention to the mode trapped in the temperature adjustment layer. It is 
easy to show with the aid of expressions (3.11) and (3.12) that in this temperature adjustment 
layer, 

uv x + vv y = -j^[BM 3 ' 2 + g- - 7 ff] + 0(1), (3.24) 

where 

B d = lim M~ z ' 2 r T {rt)drj. (3.25) 

M — *00 Jo 

An investigation of the y-momentum equation (3.19) shows that the Gortler number must be 
of order M 3 / 2 in order to enter the leading order analysis. Thus we write 

i«(i)G = G*(z)Af 3 / 2 , (3.26) 

so that for a given constant Gortler number G we compute G*(z) using 

G*(x) = (z)GM- 3 / 2 . (3.27) 
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(3.28) 


For convenience, we also define another function Q(x ) by 

Q ^ = (23)3/2’ 

so that 

uv x + vv v + i k(x)Gu 2 = ( G * - Q)M 3 / 2 + o(M 3 / 2 ). (3.29) 

With the use of this relation, we can deduce from the perturbation equations (3.18)— (3.22) 
that 

V = 0(M 3/A T), W = 0(M 3 ' a T), P = 0(M 3/2 T), U = 0{M^T), (3.30) 

and that for fixed r}, 

■k = ot*n 

where 

Mi d = 


We therefore look for asymptotic solutions of the form 

T = exp (m 3/a J P{x)dxj {T 0 (x,r]) + M~ 3/i Ti(x, 77 ) + •■•}, 

U = Mf 1 exp ^M 3 / 4 y P(x)dxj {U 0 (x,7)) + M~ 3 ^ a Ui{x,t]) -+ }, 

V = M 3 / 4 exp (m 3 ! a J* f3(x)dxj {V 0 (x,r]) + M’^V^x.t?) + •••}, (3.31) 

W = M 3 / 4 exp (m 3/4 j X p{x)dxj {W 0 (a,i?) + ilf- 3/4 Wi(x,7;) + •■•}, 

P = M 3 / 2 exp ^M 3 / 4 J P(x)dx ^ {Po(x,r;) + M~ 3 ^ A Pi(x,r]) -f }, 

where the spatial amplification has been taken care of using the WKB method and P{x) is 
the local growth rate to be determined. On substituting (3.31) into (3.18)— (3.22) and then 
equating the coefficients of like powers of M, we obtain a hierarchy of equations. To leading 
order, we find that Vq satisfies the differential equation 


9 ! V„ 2 T'SVo am „ 

° 

whilst T 0l Wo and Pq are related to Vo by 


L 2 

V2xp 2 


Q)f'V 0 , 


T 0 = - 



iaWo = 


1 dv 0 

V2 xT drj ’ 


pVlxdl o 
T 2 k 2 drj ’ 


(3.32) 


(3.33) 


and Uq does not appear in our leading order analysis. Here k d =' \[2xa is the local wavenumber. 
Equation (3.32) subject to Vo vanishing at rj = 0, oo is a Sturm- Liouville problem which has 
solutions if 


(G* - Q)T' 

P 2 


< 0 . 
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This means that 


/3 2 > 0 if G'>Q, P 2 <0 if G*<Q , 


since T 1 < 0. It then follows that neutral stability (/ 3 = 0) occurs at the position x = x n where 


G*(x n ) = Q(Xn) 


(3.34) 


at zeroth order. Therefore, in view of the definitions (3.27) and (3.28), the neutral Gortler 
number has the expansion 

2 jj 

G = — — — TgTj M 3 / 2 + higher order correction terms. (3.35) 


For the rest of this paper we shall take Gm to be the zeroth order approximation to the “critical 
inviscid Gortler number”, thus Gjy is obtained by retaining only the first term on the right 
hand side of (3.35). An important point concerning (3.35) is that the first term on the right 

—3 

hand side is independent of x n if the wall curvature varies like x s ; in the latter situation 
nonparallel effects dominate and the vortex growth rate is smaller. Thus to determine the 
higher order correction terms to the neutral Gortler number, we have to distinguish two cases, 
namely, (i). k(x) = ( 2x)~ 3 / 2 and (ii). k(x) ^ (2 x) -3 / 2 . They will be treated separately in the 
next two sections. 

On the other hand, equation (3.32) can also be interpreted as an eigenvalue problem which 
determines the growth rate f3(x) at a given value of x corresponding to any wavenumber k. 
The appropriate boundary conditions are deduced as follows. As r] — » oo, T — ► 1 and equation 
(3.32) reduces to 


d 2 V 0 

dr/ 2 


- k 2 v 0 = o, 


so that Vq ~ exp(— krj) and the asymptotic condition 


?p- + kV 0 = 0 

or) 


(3.36) 


should be imposed at “infinity”. As 77 — » 0, T — ► A 2 /r f where from (3.15b) A = 3(1 + m)/cr. 
Equation (3.32) reduces to 


d 2 V 0 8dVo _ k 2 A A (G* - Q)k 2 4A 2 
dr] 2 r] dr) rf ° y/2x/3 2 rf °’ 


which has the solution 


Hence the asymptotic condition 


v ^ hA \ 
^0 ~ exp(--^3-). 


S ^- h 4vo= 0 

dr] rj* 


(3.37) 

(3.38) 

(3.39) 


should be imposed at “zero” (which is taken to be some small value in numerical calculation). 
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The eigenvalue problem (3.32), (3.36) and (3.39) was Bolved numerically by employing a 
fourth order Runge-Kutta method. In Fig.3, we have shown the dependence of the growth rate 
on the local wavenumber. The plot clearly shows that as k — > 0, /3 2 — > 0 whilst as k -* oo, 
/ 3 2 — y constant. These features are also borne out by the asymptotic analysis given in the 
following subsection. The inviscid mode we have described above therefore has growth rate 

3 

proportional to Mi and we refer to it as the strongly unstable inviscid mode. We note that 
when G* = Q the growth rate vanishes. In this case it is necessary to look for evolution of the 
vortices on a shorter lengthscale in the streamwise direction; that problem will be addressed 
later in this paper and we shall refer to the inviscid mode in that regime as the near neutral 
inviscid mode. 


3.2 The small and large wavenumber limits of the strongly unstable inviscid 
mode 

First, we note that as k — > 0, we are approaching the scalings for the wall mode (see §6) which 
has wavenumber k ~ 0(M ~ 3 / 2 ) and which is trapped in the wall layer. Thus Gortler vortices 
are appropriately governed by (3.37). The solution (3.38) shows that vortices decay to zero in 
the thin layer rj = 0{k 1 / 3 ) near the wall. This is also verified by the obvious shift to the left of 
the first mode eigenfunctions in Fig.4 with decreasing k. It can then be deduced from (3.37) 
that p 2 has to be 0{k) in order to enter the leading order analysis. 

Next, in the large k limit, a WKB analysis of (3.32) shows that Gortler vortices will be 
trapped in an 0{k~ l l 2 ) thin layer centred at r\ — rf where fi 2 has a maximum. Thus we 
introduce a new variable ( by 

c = * l/a ( 17 -*?•), (3-40) 

expand (G* — Q)/(\/2xP 2 ) as 

^2~jp = + k 1 / 2 \i+k 1 A2 + >,> , (3.41) 

and look for solutions of the form 


VoM = F 0 °(*,0 + ~k- 1/2 v 0 \*,0 + ~k- l v 0 2 (x ,0 + • ■ • . 


(3.42) 


On substituting (3.40)-(3.42) into (3.32), equating the coefficients of like powers of k , and 
then solving the resulting set of equations, we find that to leading order, Aq is determined as 


Aq = - 


TV) 

T'(vY 


(3.43) 


At order k 1 / 2 , Ai is determined as Ai = 0 if we insist that /3 2 attains its maximum at tj = rj*. 
At order A -1 , Vq 0 is found to satisfy the parabolic-cylinder equation 


d 2 V 0 ° 

dp 


-^V o °-aF o ° = 0, 


(3.44) 
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where 


„ f^rfu * \ f 1 d 2 , 1 x l 1/4 . _ \ 2 T'(v*) f 1 d 2 . 1 .1 1/2 /0 „ c . 

i- 'T (t? )| 2 f/( 7? *) cil? ^( Ao )| a " 2f 2(7?*) { 2 ^( 7 ?*) J ' ■ 

If we impose the condition that the disturbance is confined to the fc -1 / 2 layer we must choose 
a = —(1/2) - s, where s is a non-negative integer. This condition determines the infinite 
sequence of eigenvalues 


A2 = = Aq 


1 + 23 




/-2T , (t7*)V A 0 


The eigenfunctions corresponding to the eigenvalue A 2 * are given by 

V 0 ° = = e-^He,(0, 


(3.46) 


(3.47) 


where He,(£) is a Hermite polynomial. 

With the aid of the numerical values for T given earlier in this section, we find that 

= 0.5786(1 - 0.5297^4^ + •••}, (3.48) 

(j — Q k 

and that rf = 2.3228. Both these numerical values and the concentration of vortices in an 
0{k~ x I 2 ) region are confirmed by the numerical results shown in Fig.3 and Fig.5. 

In closing this section, we note that results given here for hypersonic flows are in sharp 
contrast with related results for incompressible flows. In a recent paper, Denier, Hall and 
Seddougui (1990) have discussed the spectrum of the large Gortler number eigenvalue prob- 
lem. It was found that the inviscid Gortler vortex eigenvalue problem has an exact solution 
with the spatial growth rate increasing monotonically from zero and tending to infinity at 
large wavenumbers. In the high wavenumber limit viscous effects become important when the 
wavenumber is 0(G») and a maximum growth rate is achieved in that regime with the growth 
rate tending to zero in an 0(G*) regime as discussed by Hall (1982). At small wavenum- 
bers the growth rate tends to zero and the vortices spread out above the boundary layer. In 
fact when the wavenumber is 0(G~ r) an eigenvalue problem related to that for Tollmien- 
Schlichting waves is recovered but at even smaller wavenumbers nonparallel effects dominate 
and the problem must be solved numerically as in Hall (1983). The major difference we have 
found above for hypersonic Gortler vortices is that the growth rate tends to a finite value as 
the local wavenumber tends to infinity; we shall see later that this has a significant effect on 
the way in which viscous effects come into play at high wavenumbers. 


4 Neutral instability with k(x ) ~ (2x) 3 / 2 

We now proceed with the determination of the higher order correction terms in the neutral 
Gortler number expansion (3.35). In the case when the curvature «(x) = (2x) -3 / 2 , G*(x) = 
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Q(x ) and the 0(M 3 / 2 ) term on the right hand side of (3.29) vanishes for all x. The implication 
of this is that for this special distribution the curvature of the basic state is exactly counteracted 
by wall curvature over an 0(1) interval in x, in the more general case that is only the case over 
an asymptotically small interval. An investigation of the perturbation equations (3.18)— (3.22) 
with the aid of the equations (A1)-(A5) given in the Appendix then reveals that the neutral 
Gortler number must expand as 

G = 2BM 3 / 2 +G + o(l), (4.1) 

and that the perturbation quantities have relative orders 

P = 0(^-n T = 0(V), W = 0(V), P = 0(V), 

where G = 0(M°) is to be determined. We therefore look for the following form for the 
solutions for (3.18)- (3.22): 

u = j£U(x,v) + ‘~, V = V(x,r,) + W = W(x, V ) + •••, (4.2) 

P= -Z=P{ x > r l) + ■■■ , T = V2 x0{x,7]) + •••, 

V2x 

where the insertion of the factor \/2x is purely for convenience. 

On substituting (4.2) into (3.18)— (3.22) and then equating the coefficients of like powers 
of M , we obtain, to leading order, 


( 4 . 3 ) 


( 4 . 4 ) 


( 4 . 5 ) 


d .fldU, du - r2 ~ „ dU 

-joz{!lv + a t" - EAs - BE?*] 

V + 9|/ ( f )] f drj j 

dP rif- dV _ y2fR - ,4 d ,fidV, 
dr/~ T V + T] dr] M V + 3d7](fdi ^ 


+l[ V 2 


oT 




~ ~ ¥ T ' m + ~ 2x 


3 drj K T 
1 d n dV 

rr^ tkW) - 2X ~d: 7 
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(4.7) 




Here ib =’ \J2xa. It can be seen that (4.4)— (4.7) are independent of U and the latter is 
determined from (4.3) after ( V , W, T, P ) have been determined. Since these leading order 
perturbation equations are parabolic with respect to the variable x, they have to be solved 
by specifying the perturbation quantities at a given upstream position and then marching 
downstream. We therefore expect that neutral stability would depend crucially on what initial 
conditions we impose and where we impose them. However, before we present our numerical 
solutions of these equations, we first consider a special case: the large wavenumber limit for 
which a simple asymptotic solution is possible. 


4.1 Large wavenumber limit 

In the large wavenumber limit, the length scale over which vortices vary is small compared 
with the lengthscale over which the boundary layer grows. Then we expect that nonparallel 
effects do not come into our leading order analysis. This is indeed the case, as we show below. 

For large k, vortices are confined to a thin layer of 0(e 1(/2 ) thickness centred on 77 = 77 * 
where, 

/k, (4.8) 

and where rf is the most unstable position to be determined in the course of our calculation. 
We therefore define a new variable 0 by 

0 = e~ 1 / 2 (V~ rf). 

An investigation of equations (4.3)— (4.7) shows that when e — > 0, 

9 = 0(e 2 V), W = 0(e^ 2 V), P = 0(<T 1/2 F), 

G = 0( 1/e 4 ), A = 0 (l/e 2 ). (4.9) 

Hence we look for the following form of asymptotic solutions for (4.3)— (4.7) 

G = ^4 (00 + + €02 + e 3 ^ 2 03 + •••)> 

V = (Fo + e 1 / 2 Vi + • • -)E, W = {e^ 2 W 0 + eW 1 + • • • )E , (4.10) 

P = (€- 1 / 2 P 0 + P 1 + 9 = (e 2 9 0 + e >/2 9 1 + • • • )E , 

where 

E = e xp J (/3o(0) + e 1/2 A(0) + • • » (4.!1) 
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and where Vq,V\ etc. are functions of <f> and x. Note that E here represents the fast variation 
of the perturbation quantities along the stream wise direction whilst the dependence of Vo, Wo 
etc. on x represent the slow variation of perturbation quantities due to the nonparallel effect 
of the boundary layer growth. In effect we have described the fast variation of the disturbance 
by a WKB type of expansion in the streamwise direction. Here we are only concerned with 
neutral stability, so we set Po = P \—P 2 = 0. On substituting (4.10) into (4.3)-(4.7) and then 
equating the coefficients of like powers of e, we obtain a hierarchy of matrix equations. To 
leading order, we have 


00 ~ pJs Vot 

fio dVo 

-no = - 


iW 0 = - 


Gq = - 


1 dv 0 

T 0 d<j> ' 
2 /2§f 0 4 


T 0 d<f> ’ u oT x ’ 
where To = T(r]*),Ti = T'(r]*) and fio = fi(To). To next order, we deduce that 


(4.12) 


G 1 = 0 . 

At next order we find that Vo must satisfy the parabolic-cylinder equation 

- \? v ° - aFo = °' 


Here 


, , r in 2 Tg 

( ~* X • °- G! 6^ A ’ X --3Go[-W-, 


(4.13) 

(4.14) 

(4.15) 


If we impose the condition that the disturbance is confined to the e 1 / 2 layer we must choose 


a = -1/2 - s, s = 0,l,*", (4-16) 

The smallest G 2 corresponds to s = 0 and we then have from (4.15b) that 

G - ( 3g ° 1 u 17 n 

■ W erf J, = „. ■ ( - 17) 

The centre of vortex activity 77 * is determined by the condition that Go attains its minimum 
there: 

W = 0. (4.18) 

After solving (3.14) numerically for the basic state temperature T, we then use (4.12d) and 

(4.18) to determine 77 *, and (4.12d) and (4.17) to determine Go and G 2 . We find that 

77 * = 3.001, Go = 15.4834, G 2 = 34.3175, 


so that 

G = 15.4834fc 4 + 34.3175ft 3 + • • • . (4.19) 

Finally, we remark that the above analysis is valid as long as the local wavenumber ft = y/2xa 
is large. This means that the far downstream evolution of Gortler vortices can always be 
described by the above theory. 
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4.2 Numerical results for k 0(1) 

When the wavenumber is 0(1), the perturbation equations (4.3)-(4.7) are parabolic with 
respect to x and they have to be solved numerically subject to some initial conditions imposed 
at some upstream position x. For computational purpose it is convenient to eliminate P and 
W among (4.4)-(4.5). After some manipulation, we obtain 


where 


1 d*V 2x d 3 V 4 xT' d 2 V k 2 T 2 n dV ^ 

+ - « « JT 2x ~fa = £' i ai ’ 


T drj A fi dxdrj 2 yT dxdrj 


_fi_d 2 0 = fit' 2fif> dO T'~ 
oTdr? VT 2 aT ^ drj + T 


^oT 2 aT 

7]T' _ 1 d 
T a dry T 


dO_ 

dx 


+ I 1 + TT - 


01 = " PT 


a, = P[^ V tF- V r+ 

II fi 1 ‘dT] 


“3 




+ V T” - ^ - ^(p.'vf' + fif + fi V f") 


+ M 7? f'2 + ^ 2 + %fff' , T „ 


+ 2=n 9 ' 


an 


- { ~-W + - §£) - *’(f + ^(T"f - r°)] 


y2 


3T' 


- ^' 2 + fP(^' 2 + T' " - “r) 


+ ^7 ? rf , -i 4 f 3 |(-F), 

3f'9 3 F d 3?' d 2 V d 2 f dV 5 3 W 3 

a 5 — O ■> — a„\ rfi 2 ) a„2 XJ„2 — 


+PT 2 


f 2 d ?/ 3 dr] K f 2 dr? dr] 2 K T 2 ' drj drj 3 

\/i T+ A + T + PyT 2 ) + fi dxdr) 2x fi dx 1 
fc 2 ? 2 50 

“s = —2x7/—, 

/x ax 


a 7 = -2xTk 


r 2 5 2 0 


drjdx' 


as 


1 r fi fi J 


5 2 W 

5t/ 2 ’ 


(4.20) 


(4.21) 
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ag 


-{ 


2W yP_+2T £ 1 , ( >£ _ f „ } . Ep \ 9*. 

ft ' ft T A T ' ft J drj 


In these equations, W 3 and W are given by 
W 3 = 


ft d 2 6 + ^2ftT' 




0 T 2 drj 2 K oT 2 aT 3J dr] 


< 4 - 22 > 


W = -i— + W 3 . 
T 07/ 3 


(4.23) 


The finite difference scheme used here is similar to that used by Hall (1983), Wadey (1990) for 
incompressible flows. The reader is referred to the latter papers for a more detailed discussion. 
Our implementation of the numerical scheme is as follows. We first specify V, U and 9 at a 
given upstream position x. Then we use the finite difference scheme to march downstream and 
thus calculate the evolution of the initial perturbation with respect to the streamwise variable 
x. The position of neutral stability is defined as the place where a certain energy measure has 
zero growth rate. We use the following five energy measures to monitor the energy variation: 


Ei = U 2 dy = V2x f°° U 2 Tdrj, 

Jo Jo 

E 2 = /°°(F 2 + W 2 )dy = \PIx [ °° (V 2 + W 2 )Tdr], 
Jo Jo 

E 3 = f \plxQdy - 2x f dfdrj, £4 - Umax, 
Jo Jo 


(4.24) 

(4.25) 

(4.26) 

(4.27) 


5 ~ I dV ~ V2iJo t ? ) dv - 

We give results for several measures of the vortex strength in order to show the relatively large 
variation of growth rates associated with different flow quantities. Depending on the initial 
conditions and on which energy measure we use, the disturbance can either grow or decay 
initially. In the former case, there is only one neutral position corresponding to each pair 
(a, ( 3): the energy will reach a maximum at some downstream location and then decay to zero 
monotonically. The corresponding neutral curves only have right branches. The region on the 
left of a neutral curve is unstable and the region on the right is stable. In the latter case, if 
the Gortler number is large enough, there are always two neutral positions: the energy will 
reach a minimum at some downstream location, then grow to reach a maximum and finally 
decay to zero monotonically. The corresponding neutral curves have distinct right and left 
branches. The region above a neutral curve is unstable and other regions are stable. For the 
special curvature case discussed in the present section, the flow is always stable for x >• 1 , since 
according to the asymptotic result obtained in the previous subsection, the flow is neutrally 
stable where G ~ k A for k > 1 so that k(x) must increase at least as quickly as x 1 / 2 if the 
vortex is to be unstable for x > 1 . 
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At a neutral position, we calculate the local wavenumber a x and local Gortler number G x 
defined by 

a x = k — \^2xa, G x = /c(x)(2x) 3 / 2 (5 = G. (4.28) 

By fixing a and varying G, we obtain a series of such points (a x , G x ) and thus plot a 
neutral curve. In our calculation, a is fixed at 0.1; the step lengths along x and 77 directions 
are taken to be 0.1 and 0.05 respectively; and the lower boundary of the 0(1) temperature 
adjustment layer is taken to be 0.8 whilst infinity is approximated by 20.8. 

In Fig.6-Fig.12, we have given the results of our numerical calculations. Fig.6 shows live 
different neutral curves obtained when we use the five different energy measures (4.24)— (4.27) 
to monitor the energy growth. Fig.7 shows three different neutral curves which are obtained 
when we impose three different initial conditions at the upstream position Xq = 20, whilst Fig.8 
shows three different neutral curves which are obtained when we impose an initial condition 
at three different upstream locations. In two of these graphs, we have also plotted the two 
term asymptotic result (4.19). As we expect, although these neutral curves have distinct left 
branches, their right branches all converge to the unique large local wavenumber limit. We 
can see from Fig.8 that as the initial location of the disturbance moves towards the leading 
edge, the neutral curves move progressively up and across to the right. This is contrary to 
the corresponding results found by Hall (1983) for incompressible flows. Finally, Fig.9-Fig.12 
show the typical profiles of the four perturbation quantities as they evolve downstream, with 
the initial conditions given by (iii) in Fig.6 and Xq = 20 ,6? = 1000. In order to see how the 
the centre of vortex activity evolves downstream, we have normalized each of the perturbation 
quantities by its maximum. It is clear that as the vortices evolve downstream (and thus as the 
local wavenumber increases), they become more and more concentrated, which agrees with the 
asymptotic result found in the previous subsection that in the large local wavenumber limit, 
Gortler vortices are trapped in a thin layer of depth 0(k~ 1 / 2 ) centred at 77 = 3.001. 

5 Neutral instability with k(x) ^ (2x)~ 3 / 2 

When the wall curvature is not proportional to (2s) -3 / 2 , the 0(M 3 / 2 ) term on the right hand 
side of (3.29) only vanishes at the leading order neutral position and its effect will persist in 
the downstream development of Gortler vortices. An important consequence of such an effect 
is that non-parallel effects will be important over a larger range of wavenumbers than was 
the case for the special curvature case. Suppose we measure the order of the wavenumber by 
writing it as a = 0(M a ). Then we will show in this section that non-parallel effects continue to 
be dominant for a up to and including 1/4. For a > 1/4, non-parallel effects become negligible 
compared with viscous effects and an analytical expression can be obtained for the second 
order correction to the Gortler number expansion. This second order correction becomes of 


20 


the same order as the leading order term due to the curvature of the basic state when a = 3/8. 
In this case a higher order correction term can also be obtained. 


5.1 0(1) wavenumber regime — the near neutral inviscid mode 


In the 0(1) wavenumber regime, it is convenient to determine the stability properties by 
considering the evolution of Gortler vortices in the neighbourhood of the leading order neutral 
position x n given by (3.34). Thus we shall fix the Gortler number as given by 


G = 


2 B 

(2 x n ) 3 / 2 K(x n ) 


M 3/2 


(5.1) 


and determine the second order correction say x n to the neutral position x n so that Gortler 
vortices with G given by (5.1) are neutrally stable at location x n + x n . Replacing x n by x n - x n 
in (5.1) then gives the appropriate expansion of the Gortler number for vortices neutrally stable 
at x — — Xfi . 

It can be shown that in the neighbourhood of x n , the second term in the expansion of 
k(x)G/ 2 will force a growth rate of order M 1 / 2 . Hence we shall consider the evolution of 
Gortler vortices in an 0(M -1 / 2 ) neighbourhood of x n by defining a new variable X by 


X = (x- x n )M 1 < 2 , 


(5.2) 


and look for asymptotic solutions of the form 

T = T 0 (X , t?) + • • ■ , V = M^ 2 V 0 (X, t?) + ■ • • , 

W = M 1 / 2 W 0 (X, t?) + • • • , P = MP 0 (X, 77) H . 

Equation (3.29) becomes 


where 


uv x + vv y + -k(x)Gu 2 = EXM + o(M), 


n def. d(G*-Q), 

E ~ dx U= *"- 


(5.3) 


(5.4) 


(5.5) 


Note that it is this term that gives rise to a local growth rate of order M 1 / 2 . On substituting 

(5.3a^d) into the perturbation equations (3.18)-(3.22), dropping higher order terms and then 

eliminating Vo, Wo and P 0 from the resulting equations in favour of T 0 , we find that Vo, Wo 

and Pq are related to To by __ 

s/2 x n T dT 0 


Vo = —■ 


ikW 0 = - 


t' ax’ 

1 dV 0 


Pq — 


T dr) ’ 

s/2x^ a 2 Vb 
fc 2 T 2 dXdrj 2 ’ 


(5.6) 

(5.7) 

(5.8) 
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and that Tq satisfies the equation 


a 2 f a 2 r 0 

dX 2 \ dr] 2 


2T" OTq 
T' drj 


- [& a r : ! - Tf'(± - ^)']r 0 } - k 2 


EXT' 

V^ii 


T o = 0, 


(5.9) 


where k d =‘ y/2x n a. We can interpret (5.9) as the turning point equation associated with the 
breakdown of the WKB structure in x of the expansions (3.31), indeed the evolution equation 
(3.32) is retrieved from (5.9) by taking X to be large. Since this breakdown is associated 
with a simple zero in x of the right hand side of (3.32) we expect that the local behaviour of 
the disturbance should now be expressed in terms of Airy functions. Equation (5.9) admits 
separable solutions of the form 

T 0 (X,r,) = i>(X)iP(r,), (5.10) 


with <j> and ip satisfying 

4>"{x) - U)X4>(X) = 0, (5.11) 

1 'pff FT* 1 

r(v) - ^'(v) - [ k 2 f 2 - TF(± - - k 2 • = o, (5.12) 

where the separation constant u is to be determined by solving the eigenvalue problem (5.12) 
subject to appropriate boundary conditions. By a simple substitution z = Xw 1 / 3 , equation 
(5.11) reduces to the standard form of Airy’s equation W"(z) - zW{z) = 0 which has two 
independent solutions Ai(z) and Bi(z), so the solution of (5.11) is given by 


4>{X) = aAi(u>^ 3 X) + bBi{u l < 3 X), 


(5.13) 


where a and b are two constants to be determined by initial conditions. 

To solve (5.12), we first note that in the large wavenumber limit, equation (5.12) takes 
the same form as equation (3.32). Therefore, the solution of (5.12) can be written in terms of 
Hermite polynomials as 

i’iv) = i>, = • He,{£) (5.14) 


and from (3.41) the eigenvalue u> expands as 







(5.15) 


where Ao, A 2 3 and £ are defined in turn by (3.43), (3.46) and (3.45a). 

In the 0(1) wavenumber regime, equation (5.12) has to be solved by a numerical integration, 
and in general an infinite number of eigenvalues u a (s = 0, 1, • • •) and eigenfunctions ip, can be 
obtained. Then the general solution of (5.9) can be written as 


T 0 = Y^ { a„Ai{u) a X ) + b,Bi(u 3 X)}ip a (r]), 

s = 0 


(5.16) 
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where a , and b a are constants to be fixed by initial conditions at X = 0. From (5.6), Vo is 
given by 

Vo = £ {a.Ai'{u,X) + b 3 Bi'(u a X)}cj t Tp,(r,). (5.17) 

1 4 = 0 

It is clear that once To(X,rj ) and Vo(X,ri) are specified at X = 0, the coefficients (a 4 ,& 4 ) and 
hence the evolutionary behaviour of the perturbation field (Po, Wo,To,Po) will be completely 
determined. 

The correction term to the neutral position can be defined as the position where a certain 
energy measure has zero growth rate. It is obvious that such a position would depend upon 
what initial conditions we impose at X = 0 and what energy measure is employed to monitor 
the energy growth. In principle then it is an easy matter to determine the local neutral position 
associated with any initial perturbation, we note however that before growth of the vortices 
occurs they will have an oscillatory behaviour in X since both Airy functions axe oscillatory 
on the negative real axis. Clearly this occurs because the boundary between instability and 
stability is controlled by inviscid effects in this regime, there is no counterpart to this result in 
the behaviour of Gortler vortices or for that matter Tollmien-Schlichting waves in incompress- 
ible flows. We further note that appropriate forms for the initial conditions can be obtained 
from the receptivity problems associated with wall roughness or free stream disturbances, see 
Denier, Hall, and Seddougui (1990) and Hall(1990). We merely note in passing here that it is 
reasonable to expect that the type of mode discussed above is more likely to be stimulated by 
free-stream disturbances since the effect of wall roughness is diminished by the wall layer over 
which the wall roughness must diffuse before reaching the unstable adjustment layer. 

This is certainly typical of the evolution of Gortler vortices in growing boundary layers. In 
the present problem, non-parallel effects dominate in the evolution of Gortler vortices mainly 
through the 0(M 3 / 2 ) curvature of the basic state. As we increase the wavenumber, viscous 
effects will gradually come into play in the evolution of Gortler vortices and nonparallel effects 
will become less important. In the following subsection we consider wavenumbers of order 
M 1 / 4 . This is the maximum order at which nonparallel effects are dominant. We shall show 
that when the wavenumber is increased further above this order, nonparallel effects become 
negligible. 

5.2 The 0(M X / 4 ) wavenumber regime — the nonparallel viscous mode 

When the wavenumber reaches the order M 1 / 4 , the streamwise lengthscale is still 0(M -1 / 2 ) 
(implying that the local growth rate is 0(M 1 / 2 )), but the vertical lengthscale becomes of 
0(M -1 / 8 ) (as we expect that vortices would be trapped in an 0(M -1 / 8 ) thin layer). We 
therefore define two new independent variables X and £ by 

X = (x - x n )M'l\ £ = (*?- (5.18) 
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where rf is the center of vortex activity. For convenience, we define a small parameter e by 




(5.19) 


An analysis of the perturbation equations (3.18)— (3.22) shows that 

V = 0(e~ 4 T), W = 0(e~ 3 T), P = 0(e~ s T). 


(5.20) 


We now assume that T = 0(1) and scale ( V,W,P ) by (e 4 ,e 3 ,e 5 ). Then by neglecting 
terms of relative order higher than e 2 in the perturbation equations (3.18)-(3.22), we obtain 

4 


dX + M T ' 3-(2x n )d{ K T d{ 


d ,fidV 

'•'t ac } 


+ e 2 


1 dP 

V 2^ dZ 


dW 


— € 


1 dW - 4 - - 

- TP + -pTW + 


3x/S^ d( 
p dV 


= 0 , 


(5.21) 


d 2 dX 


+ €■ 


pT' 


3\/2x n d£ \/2x n 


:V 


— e 


d (£°W) = o 

at ) u > 


(2x n )a 2 dCT dZ 


W = - 

V T dT 


1 dV 1 dT 
y/2^T dZ +€ ' TdX + e ' 


—V 

s/2T n T 2 ’ 


a 2 pT 2 


T + e 2 


d m( H9T_\ 

at A 


(5.22) 

(5.23) 

(5.24) 


s/2x, ~ T'dX oT' {2x n )oT' dZ K T dZ 

where W d = iaW , a d =’ Af -1 / 4 a and where we have used the same notation to denote the scaled 
perturbation quantities. Elimination of V, W, and P in favour of T among these equations then 
gives 

d 2 T 1 . FIT 1 . F.YT \ 

r T 


S ’ r i (. , ‘VA7T 91 I >T EX 

+ (1 + -)a + -a T„ T + ^=jT 5 - 


= 0 . 


(5.25) 


where ?o = T{t]*),T\ = T'(r]*),po = p(To ) and where 77 * is chosen at higher order such that the 
vertical structure of T can be expressed in terms of parabolic cylinder functions which vanish 
at Tj = ± 00 . The solution of the above equation is easily found to be expressible as the product 
of an Airy function multiplied bye 2 . As before the Airy function grows or decays 

exponentially for large X and is oscilatory on the negative real axis. However the presence of 
the exponential factor now leads to a crucial change in the nature of the streamwise evolution 
of the disturbance. We refer to the fact that the exponential factor, induced by viscous effects, 
now means that in the stable regime the disturbance decays exponentially rather than oscilating 
as was the case previously. This result is consistent with the usual result of stability theory 
that inviscid disturbances change from being oscilatory to being exponential in character when 
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instability occurs; viscous instabilities on the other hand are exponential in nature either side 
of the stability boundary. The local neutral position can only be obtained by specifying an 
initial disturbance and finding where the solution of the evolution equation begins to grow. 

If we let a<l, then to leading order equation (5.25) reduces to 

d 2 T EXTi 


dX 2 


+ 




T = 0 


(5.26) 


which matches with the asymptotic form of (5.9) when k > 1. 

On the other hand, if we let a > 1, then to leading order (5.25) reduces to 


laVT 'A T = 0, 

Vv^Ti 2 <r J 


(5.27) 


which shows that the second order correction term to the leading order neutral position becomes 
independent of nonparallel effects and is given by 


X = X n = - 


- 


crETi 


r ,2rpi~i 

HqIqO. , 


(5.28) 


On replacing x n by x n — M~ l ! 2 X n in (5.1) and expanding the resulting expression in the 
neighbourhood of x n , we obtain the following expression for the Gortler number for Gortler 
vortices neutrally stable at location x = x n : 

2 BM*! 2 


G = 


+ o 4 <7o + 


where 


k(x„)(2x„) 3 / 2 

2\/2x n p, 2 TQ 


9o = -- 


PKqTi 


(5.29) 


(5.30) 


The expression (5.29) is valid for wavenumbers of order M“,l/4 < a < 3/8. When the 
wavenumber reaches the order M 3 / 8 , the second order correction term in (5.29) is as large as 
the first term and a more accurate asymptotic expression can be obtained. That situation will 
be discussed in the following subsection. We have seen above that in the 0(M« ) wavenumber 
regime viscous effects come into play and modify the evolution of the near neutral inviscid 
mode, for that reason we refer to the mode in this case as the nonparallel viscous mode. We 
further note that the strongly unstable inviscid mode connects directly with the parallel viscous 
mode discussed below. Thus we have shown above that the nonparallel viscous mode connects 
with the near neutral inviscid mode in the vicinity of the right hand branch of the neutral 
curve. Hence the initial stages in the evolution of the right hand branch of the neutral curve 
are governed by an interplay between viscous and inviscid effects. We further note that, in 
view of the limiting form (5.29) valid for large a, as the wavenumber increases the neutral 
Gortler number will also increase. 
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5.3 The O(MS) wavenumber regime — the parallel viscous mode 


When the wavenumber becomes of order M 3 / a , viscous effects are of the same order as the 
centrifugal acceleration of the basic state in the determination of the Gortler number, and the 
leading order inviscid result (3.35) is to be modified. We assume that to leading order the 
Gortler number now expands as 


G = 


2BM 3 ' 2 , 4 

*(*n)( 2x n )3/2 + “ 9 ° * 


(5.31) 


Here the first term is due to the curvature of the basic state and the second term is due to 
viscous effects and is to be determined. 

For convenience, we introduce a small parameter e and an 0(1) constant N by 


e d ^~, N d ^M 3 ! 2 e\ 


so that (5.31) can be written as 


__2 BN__ 1 so 
K (x n )(2x n f! 2 ' e^‘ 


(5.32) 


(5.33) 


To determine the higher order correction terms to the Gortler number expansion, we shall first 
fix the Gortler number as given by (5.33) and consider the evolution of Gortler vortices in the 
neighbourhood of the leading order neutral position x„ defined by (5.31), aiming at finding the 
second order correction say ex n to the neutral position. As we have remarked at the beginning 
of the first subsection, replacing x n by x n — ex n in (5.33) would give the appropriate expansion 
of the Gortler number for vortices neutrally stable at x = x n . 

The vortices under consideration vary on small lengthscales in both x and Tj directions. 
In the stream wise direction, their growth rate can be shown to be 0( 1/e) so that they evolve 
on an 0(e) lengthscale. In the rj direction, they are confined to an 0(e 1 ^ 2 ) thin viscous layer 
because of their small wavelength character. We therefore define two new variables X and ( 
by 


X = 


X — X r 


sV-V 

— i /■> > 


e • ’ £ i/a 

where 77 * is the centre of vortex activity and is to be determined. 
We now look for asymptotic solutions of the form 


(5.34) 


T = Oq(X, 0 + e 1 / 2 0 1 (X, C) + e0 2 (X, () + • • • , 


V = e~ 2 [V 0 (X, 0 + e l ! 2 Vi(X, Q + eV 2 (X, () + ■■ •], 

W = e- 3 / 2 [W 0 (X,0 + e^W^X, () + eW 2 (X,() + •••], (5.35) 

P = e~ 5 / 2 [Po(X, 0 + e l / 2 Pi(X, 0 + eP 2 (X, () + ■■ •]. 
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Here the relative orders of the perturbation quantities are deduced from the perturbation equa- 
tions (3.18)-(3.22). On inserting these expansions into the perturbation equations (3.18)-(3.22), 
expanding all coefficients there about x = x n and 77 = 77 *, and then equating the coefficients 
of like powers of e, we obtain a hierarchy of equations. To leading order, the Gortler number 
go in (5.31) is determined as a solvability condition for (Vo,0o) and i 8 given by 

_ 2\/2z^qTq _ 

90 (JKqTi 

whilst 0 O , Wq and P 0 are related to Vq by 



60 — - 




2/xoT 0 2 i„ 


Vo, iW 0 =- 


dv 0 


\/2x n T 0 d( ’ 


Po = - 


f-io dV 0 

%~d( 


(5.37) 


where To = T( 77 *), f\ = T'(tj*), p, 0 = p.(To) «o = «(xn). Note that (5.36) is of the same 
form as (5.30), as we would expect. 

To next order, we obtain three expressions similar to (5.37) for $ 1 , W\ and Pi in terms of 
V\ and Vo and the condition that 

^|„ = „. = 0, (5.38) 

which implies that 77 * is where go attains its minimum. 

If we carry on one order higher, we find from a solvability condition for (V 2 , # 2 ) that Vo 
must satisfy the evolution equation 


d 2 V 0 

d( 2 


2(1 + c)Tpx n dVp 
3po dX 


- a( 2 V 0 + bXV 0 = 0, 


where 


a = 


T 0 2 x n d 2 g 0 

3g 0 


>0, 


7 _ 2j n Tp r 2 BN . K\ 3 . Ki 11 

3 lpo«o(2x„) 3 ' 1 ' 2 ko 2x„ kq 2 x„/’ 

and where k\ = n'(x n ). The solutions of (5.39) which satisfy the conditions Vo 
can be written as 


(5.39) 


(5.40) 
0 as Id — > 00 


Vb = Vb m (X,C) = exp 


3fiob 

4(1 + a)ToX r 




(da) 1 /^), m = 0,l,2,..., (5.41) 

where U is a parabolic cylinder function. The neutral position x n can be taken to be the point 
where dVo/dX = 0, so that the 7 nth mode is neutrally stable at 


The most unstable mode corresponds to m = 0. We therefore have 


x 


n 





(5.43) 


With the expression for x n determined, we can now replace x n by x n — ex n in (5.33) and then 
expand the two terms on the right hand side up to and including the 0( 1/e 3 ) term, hence 
obtaining the expansion 


_ 2BM 3 / 2 g 0 1 1 l3g 0 d 2 go 

K (x n )(2x n )3/2 + + C 3 • 2Tq dr,* 2 

for Gortler vortices which are neutrally stable at position x n . 

As is to be expected, the Gortler number expansion (5.44) agrees in the special case 
k = (2x) -3 / 2 with the combination of (4.1), (4.10a), (4.12d) and (4.17) (note that the small 
parameter e there corresponds to e/\/2x n here). This means that in the large wavenumber 
limit, the relation (5.44) is a universal expression for the neutral Gortler number, which is 
valid for all wall curvatures. 

Finally in this subsection we stress that a more unstable version of the parallel viscous mode 
can be obtained by taking go bigger than the value given by (5.36), in that case we must allow 
for a growth rate of order M ? and then (5.36) is replaced by an equation to determine that 
growth rate. This structure then enables a direct connection between the strongly unstable 
inviscid mode and the parallel inviscid mode at relatively high Gortler numbers with G—Gjj ~ 
Aft. The analysis for this more unstable regime is essentially identical to that given by Denier, 
Hall, and Seddougui (1990) in the incompressible case. 



6 The wall mode 

It has been established in section 3 that as the wavenumber becomes large, Gortler vortices 
become increasingly trapped in the 0(1) temperature adjustment layer. Thus the preceding 
three sections are devoted to Gortler vortices which have wavelength of 0(1) or smaller and 
which are trapped in the temperature adjustment layer. Clearly it is possible for vortices 
of wavelength smaller than the thickness of the transition layer to be excited, far enough 
downstream the local wavenumber will become comparable to the adjustment layer thickness 
and the previous analysis will apply. However before this occurs the vortices must be described 
by an analysis which takes account of the fact that they are of wavelength much larger than 
the adjustment layer thickness, we shall now address that situation. In fact it can be deduced 
from the definition (2.11) and (3.1c) that the variation dy of the physical variable y and the 
variation of the similarity variable dr, satisfy 

dy = V2xTdr,. (6.1) 
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The wall layer which corresponds to rj = 0(M~ 1 ^ 2 ) with T = 0(M 2 ) is therefore actually of 
0(M 3 ! 2 ) thickness in terms of the physical variable y, whilst the temperature adjustment layer 
is still of 0(1) thickness. Thus a natural scale for larger wavelength vortices is provided by the 
thickness of the wall layer, the appropriate size of the Gortler number is found by rescaling 
the vortex wavelength and velocity field by the scales relevant to the wall layer. Such Gortler 
vortices are referred to as the wall mode and are studied in the present section to complete 
our stability analysis. 

Since in the large Mach number limit the boundary layer thickens by 0(M 3 / 2 ), we should 
rescale ( y,z ) by a factor M 3 / 2 and the corresponding velocity components likewise. This 
effectively replaces all ,l R~ 1 / 2n s by u R~ l l 2 M 3 ! 2 ” . It is therefore appropriate to rescale the 
Gortler number G and the wavenumber a by defining 



G = M 3 / 2 /c(x)(2 x) 3 / 2 G, k = V2 xM 3 / 2 a. 

(6.2) 

In the wall layer, the basic state is from (2.12) and (3.5) given by 



M 3 ! 2 , 

u = F'(Y), v = — [ TF(Y) + F'(Y)n(r)], 

V 2z 

(6.3) 

where 

y 



fi(Y) d = / TdY. 

(6.4) 


'o 


The various partial derivatives of u,v and fi which appear in the perturbation equations 
(3.18)-(3.22) have to be computed before we can deduce the relative orders of the pertur- 
bation quantities. Such expressions axe given in the Appendix to this paper. With the aid of 
these expressions, we can show from (3.18)— (3.22) that the relative scalings of the velocity, 
pressure and temperature disturbance fields axe given by 

V = 0{M 3 ' 2 U), W = 0(M 3/2 U), T = 0(M 2 U ), P = 0(MU). (6.5) 


We therefore look for solutions of the form 


U = U(x,Y ) + ••., V = M 3 ' 2 V(x,Y) + ---, 

W = M 3/2 W(x,Y ) + •••, P = M 2 P{x,Y) + ---, (6.6) 

T = M 2 6(x,Y) -1 . 


On substituting (6.6) into the perturbation equations (3.18)— (3.22) and then equating the 
coefficients of like powers of M , we obtain to leading order the following set of partial differential 
equations which govern the evolution of Gortler vortices in the wall layer: 
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( 6 . 11 ) 
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Obviously, this set of partial differential equations have to be solved numerically to deter- 
mine the evolution properties of Gortler vortices in the wall layer. We further note that the 

downstream velocity component of the perturbation now does not decouple from the other 
disturbance quantities. For the purpose of numerical calculation, it is convenient to eliminate 
the pressure perturbation P and W from the above five equations. After some manipulation, 
we obtain 
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The above equations are to be solved subject to the following boundary conditions: 

At Y = 0, 

80 

•^p = 0, if the wall is thermally insulated (6.15) 


As y — > oo, 


0 = 0 , if the wall is under cooling. 


U -> 0 , V - 0 , w - 0 , 0 -» 0 . 


(6.16) 


The precise large V decay behaviour of the perturbation quantities can be deduced from 
the perturbation equations (6.12)— (6.14), but the derivation is tedious. However, the far 
downstream limit of such decay behaviour can be deduced very easily from the condition that 
the solutions of the wall mode under consideration in the double limit k — » oo and Y — > oo 
should match with the inviscid solutions (3.32) and (3.33) when k — » 0 and 77 — ► 0 there. In the 
limit k —* 0, equation (3.32) reduces to equation (3.37) and the latter has solutions confined 
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to a thin layer 77 = O^k 1 ^ 3 ) near the wall. It can be shown that in the further limit 77 — > 0, the 
solution of (3.37) has 

V 0 % (6.17) 

n 

and from (3.33), 

T 0 _ W 0 ^ (6.18) 

where ci,cj and C 3 are all functions of x. By the above-mentioned matching condition, the 
relations (6.17) and (6.18) with 77 replaced by Y are also the far downstream, large Y decay 
behaviour of the perturbation quantities in the present problem. 

For any initial disturbance located at a given downstream position, we can determine the 
position at which the disturbance becomes neutrally stable by integrating these equations using 
the same marching procedure as the one used in solving the system of equations (4.3)-(4.7). 
As an illustrative example, we impose the following initial disturbance at the location x = 50: 

(50, Y) = Y 5 e - r ' , V(50, Y) = #(50, Y ) = 0. (6.19) 


The wall curvature is taken to be y/2x and the wall is assumed to be thermally insulated. For 
a given wavenumber and a given Gortler number, we can march downstream until we reach 
the neutral position where a certain energy measure has zero growth rate. We then calculate 
the local wavenumber k and the local Gortler number G at the neutral position. By fixing 
the wavenumber o(= k( 2x)~ 1 / 2 ) at 1 and varying the Gortler number M~ 3 / 2 G from 0.0001 to 
0.03, we obtain a series of neutral points ( k , G). Fig.13 shows the neutral curves corresponding 
to the following three energy measures: 


Ei 


(°° U 2 dy = V2 xM 3/2 r U 2 TdY , 
Jo Jo 


( 6 . 20 ) 


E 2 = [ °°{U 2 +V 2 + W 2 )dy = y/2 xM 3 ' 2 [°°(U 2 + V 2 + W 2 )TdY, 

Jo Jo 

„ r°° f dU. 2j M- 3 / 2 1 .du. 2 

So V w/o f ( 9Y ) dY - 


( 6 . 21 ) 


( 6 . 22 ) 


Fig.14-Fig.17 shows the downstream evolution of Gortler vortices corresponding to the above 
conditions with M~ 3 ! 2 G = 0.001. We observe that all neutral curves decrease monotonically 
with respect to the local wavenumber and that Gortler vortices become increasingly more 
and more shifted to the right (i.e. towards the temperature adjustment layer) as they evolve 
downstream. This is certainly to be expected since in the large local wavenumber limit the 
wall mode has to match with the mode trapped in the temperature adjustment layer. By 
talcing the large wavenumber limit of the wall mode equations we can show that the neutral 
curve should tend to the limit (5 = 2B = 1.2543. In order to realize this limit numerically we 
have to carry out our calculation at very large wavenumbers. This presents some numerical 
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difficulties, because on one hand, if we fix the Gortler number and vary the wavenumber, 
the initial conditions soon become incompatible with the differential equations since for large 
wavenumbers Gortler vortices have to be trapped in a region away from the wall; on the 
other hand, if we fix the wavenumber and vary the Gortler number, a large local wavenumber 
corresponds to a large downstream neutral position and there contamination from the “finite 
infinity” becomes important because of the algebraic decay behaviour (implying that we have 
to choose a larger infinity). Finally, we remark that the ^-derivative of 6 in the expressions for 
W 3 and W could be eliminated with the aid of equation (6.13) to give a different formulation 
and thus to provide a check on our numerical scheme. We have done so and have obtained 
identical results. 

We conclude this section by stating the most important results of our investigation of the 
wall mode. We have shown that the wall layer can support a disturbance trapped in the wall 
layer with wavelength comparable with the wall layer thickness. This mode is dominated by 
nonparallel effects and has a neutral Gortler number which is a monotonic decreasing function 
of the vortex wavenumber. In the limit of high vortex wavenumbers the mode takes on a struc- 
ture essentially identical to that found for the small wavenumber limit of the inviscid modes of 
wavelength comparable with the adjustment layer thickness. Moreover in this limit the vortex 
has a neutral Gortler number which approaches from above the zeroth order approximation to 
the neutral Gortler number of the wall modes. 

7 Real gas effects 

In our previous discussions, we have assumed that the fluid under consideration is an ideal gas 
undergoing no dissociation. Our asymptotic analysis based on the large Mach number limit 
has yielded some revealing results about the stability properties of hypersonic boundary layers. 
However, in the large Mach number limit, we would expect that the wall temperature should 
be well above the temperature at which dissociation takes place. Take a boundary layer over 
a thermally insulated wall as an example. The temperature at the wall is given by 

T v , = ±{ 1 -l)M 2 T 00 

when the Prandtl number is unity. At a standard altitude of 53 km, the air temperature Too 
is 283K. If we take 7 = 1.4, M = 25, then T w = 35,375 K. Since at the far lower temperature 
of 2500 K, the oxygen molecules in the air have already begun to dissociate, it is clear that an 
investigation which takes gas dissociation into account is vitally important! 

A complete theory on real gas effects should at least incorporate the following important 
chemical reactions: 

0 2 - 20, N 2 - 2 N, 
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N + O^NO, N + O^NO + + e~. 


However, in order to expose the most important features and at the same time keep the algebra 
to a minimum, we begin our investigation by eliminating the less essential complications which 
arise from the detailed composition of air and confining our attention to a pure dissociating 
diatomic gas. The dissociation process is therefore denoted by 


A 2 ^ A + A. 


(7.1) 


In our following discussion, 0 2 and N 2 will be chosen for numerical illustration. Furthermore, 
we assume that the fluid under consideration is an ideal dissociating gas which satisfies the 
relation 

0/2 _ Pd _ T_ -T d /T 

1 - a 2 p Td 

(see Becker (1968), p. 36, or Lighthill (1957), p.6). Here a, p and T have the same meanings 
as defined in §2, whilst p j and Tj are respectively the characteristic pressure and temperature 
for dissociation. 

For the gas mixture of A and A 2 , Dalton’s law gives 

n\$lT 7i 2 SftT 
P = Pi +P2 = y ■ + —y— I 




where pi is the pressure which the ith component would exert individually if alone in volume 
V at temperature T, 7ij are mole numbers and 3? is the universal gas constant. Here we use 
subscript 1 to signify component gas A 2 and subscript 2 for A. Assume that the weight of a 
mole of gas A is m, the molar weight of gas A 2 is then 2m. We therefore have 

Til • 2m + n 2 • m 

p = v ’ 

ti 2 * m n 2 

a = = . 

Tii • 2m + 7i 2 • m 27ii + n 2 

With the use of (7.4) and (7.5), we can rewrite the equation of state (7.3) as 




p = (1 + a)9tpT, 


(7.6) 


where SR = SR/(2m) is a gas constant. 

On inserting (7.6) into (7.2), we obtain 

JL- = £l e -T 4 /T (7.7) 

1 — a p 

where pd = Pj/(KT e j) is the characteristic density for dissociation, pd, Pd and Td are in 
general functions of T, but their variations over a large temperature range are very slight. 
Their typical values are given in Table 2 (taken from Lighthill (1957)). 
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Table 2 


T d (°K) P t f(atoms) p d (gm/cm 3 ) 

Oxygen 59000 2.3 x 10 7 150 

Nitrogen 113000 4.1 x 10 7 130 

It should be noted that for atmospheric values of p , p d /p is at least 10 5 . Thus although 
p d , T d and p d are called characteristic quantities, they are not typical of the actual values 
of p, T and p at all. For p d / p = 10 5 , (7.7) shows that a is 0.05 (5% dissociation) when 
T/Tj = 0.057, and a is 0.95 (95% dissociation) when T/Tj = 0.116. For densities typical of 
the upper atmosphere, with (say) p d / p = 10 7 , these values of T /T d would be reduced to 0.045 
and 0.076 respectively. For a fixed value of p d /p = 4 X 10 s and Tj/Too = 4 x 10 3 , the variation 
of a with respect to T/T^ is plotted in Fig.18. 

Relation (7.7) was first obtained by Lighthill (1957) from quantum mechanics, and there- 
fore the ideal dissociating gas discussed here is also called the Lighthill gas. Using a purely 
mathematical argument, Becker (1968) has shown that the Lighthill gas is a special case of a 
more general class of gases. An important property of Lighthill gases is that A 2 and A have 
the same specific heats at constant volume, that is, 


Cyl C v 2- 


(7.8) 


7.1 Constitutive properties of a dissociating gas 

In order to see the complications which arise from gas dissociation, let us first recall that an 
ideal gas has the following properties which have greatly simplified our previous analysis: 

(i) The equation of state has the simple form P = $tpT as compared with (7.6); 

(ii) Specific heats c p and c v can be taken to be constant; 

(iii) The coefficient of thermal conduction k is proportional to the shear viscosity p and the 
Prandtl number a = pcp/k is usually taken to be constant; 

(iv) The viscosity can be taken to be related to the temperature by Chapman’s law, or more 
accurately, by Sutherland’s law. 

When part of an ideal gas has been dissociated, the gas becomes a mixture of two component 
gases and all of the above properties are changed. We have shown above that the equation of 
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state is modified to the form given by (7.6). In this subsection, we derive the corresponding 
expressions for /x, k and c p . Our derivation is based on work discussed in “ The mathematical 
theory of non-uniform gase s” by Chapman and Cowling (1970). 


The viscosity fi 

In the kinetic theory of gases, different expressions for the transport properties /x and k have 
been found depending on what model is used for the interaction of the gas molecules. For 
example, for a simple gas which consists of smooth rigid elastic spherical molecules, /x is 
proportional to VT, whilst if the gas is taken to consist of smooth rigid elastic spherical 
molecules each of which is surrounded by a field of attractive force, /x is given 




(7.9) 


where c is the diameter of the molecule, 5 the potential energy of the mutual attraction of two 
molecules when in contact, m the molecular weight. Relation (7.9) is known as Sutherland’s 
viscosity law. 

For a binary mixture of gases, Wilke’s law gives an approximate expression for the viscosity 
of the mixture: 


/x = 


Xi 


£L -L 
mi Mia 


+ 


X 2 


£l | £? 1 

M 31 ' M2 


(7.10) 


where /Xi and /x 2 are the viscosities of the two component gases; x\ and x 2 denote the propor- 
tions by volume of the two gases in the mixture. Since by Avogadro’s law equal volumes of 
different gases contain an equal number of molecules or moles, we have 


Til j. — u 112 t —t 1 1 \ 

*1 = ; = — ) *2 = ; = 7 - — • (7.11) 

n i + 7*2 1 + a ni + n 2 1 + a 

In (7.10), /x 12 and /X 21 are the mutual viscosities of the two component gases and are given by 

[ftr(mi + m 2 )-| 1/2 


a 


n 2 


2a 


Ml2 = 


5mi 

16c 2 


12 


27rmim 2 


/(l + -jr) i A* 21 = ~Mi 2 = (7.12) 


where ci 2 can be taken to be (c\ + c 2 )/ 2 , whilst the expression for S\ 2 has to be found empir- 
ically. Lindsay and Bromley (1950) suggested that Si 2 = \ZSiS 2 . Here Ci and Si (i=l,2) are 
constants appearing in (7.9), associated with /x,\ 

With the aid of (7.9) and (7.12), we can now write down the appropriate expressions for 
the viscosities of the component gases. These are 


fXi = 


Here 


AiT 3/2 
T + Si’ 


Ai = 


i = 1,2, /X 12 = 


A 3 T 3 / 2 


, /x 21 = /X 12 / 2 , 



A 3 = 



(7.13) 


(7.14) 


where we have made use of the fact that m\ = 2 m 2 . 
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Coefficient of thermal conduction k 


For a pure gas, the coefficient of thermal conduction is related to the viscosity fi by the simple 
formula 

k = fc v fi, (7.15) 

where / is a constant and 7 // is usually defined as the Prandtl number. 

For a binary gas mixture, Wassiljewa’s formula gives 

X\ x 2 


k = 


(7.16) 


+ ^ fr + r 

where ki and k 2 are the thermal conductivities of the two component gases and k\ 2 , and k 2 i 
are the mutual thermal conductivities. 

It is suggested in Chapman and Cowling (1970, p.256) to put 


ki _ Du 
1 ^12 j-. , 

* 12 7^12 


k 2 D 22 

I — = <* 21 - 
«21 


D\ 2 ’ 


(7.17) 


where the coefficients a * 2 and a 2 i are regarded as functions of mi/m 2 , which are determined 
semi-empirically, and for Sutherland’s model the diffusion coefficients Du, Du and D 22 are 
given by 

6 6 6 

Du = ■= — Pi, Du = - — fiu, D 22 = - — fi 2 , (7.18) 

opi 5pi op 2 

where p, is the density of the ith component gas, when pure, at the pressure and temperature 
of the actual mixture. 

If we take au = <*21 = 1 for simplicity and further make use of the above relations, relation 
(7.16) then becomes 

Pi , (f2/fl)p2 


~ Cvfl | 1 


+ 3L . i£l_ 
' *1 ma 


+ 


1 4 . 2x1 . -H 2 - 
' *a ' Mia 


(7-19) 


def. 


where we have used the fact that for an ideal dissociating gas, c„i = c „2 =’ c„; whilst fi and 
f 2 are the constants appearing in (7.15) corresponding to gases A 2 and A respectively. 

For easy comparison, we rewrite (7.10) here, making use of the relation (7.12), 

Pi , P2 


P = 


1 4 - *1 . Jil. 
si Mia 


+ 


1 2SX. . JiZ. ’ 


(7.20) 


Comparing (7.19) with (7.20) shows that we can almost write k = c v fifi for the gas mixture, 
only if f 2 — f\. It is therefore desirable to investigate the values of fi (for diatomic gases) and 
f 2 (for monatomic gases). 

For all smooth spherically symmetrical molecules, it has been shown that taking f 2 = 2.5 
provides a very good approximation. However, for diatomic and polyatomic gases, a variety of 
expressions have been suggested for /. One of these is Eucken’s formula: 


/ = ^(97 - 5). 


(7.21) 
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This formula is valid under the assumption that the transport of momentum and translational 
energy is unaffected by the internal molecular motions, and that internal energy is transported 
at the same rate as momentum. 

For diatomic gases, there are 3 translational degrees of freedom, two rotational degrees of 
freedom and on the average one vibrational degree of freedom (may vary between zero and 
two depending on whether vibration has been fully excited or not). Thus the internal energy 
is given by 

1 

Cl = 2 

(Recalling that 2m is the molecular weight of A). Consequently, the specific heats are 


(i)P +2+1 )=^' 


(7.22) 


Cvl — 


de\ 

~dT 


33? 

2m’ 


3 ? 

Cpl “ Cvl + 2m‘ 


We then have 


7 = ~ = t = !-33, / = 7 = 1-75. 


7 

Cyi 3 ' 4 

At lower temperatures when the vibrational mode is not excited, 


5,R T. 5. SR . 7.3?. 

61 ~ 2 2m 01,1 ~ 2 2m 1 Cpl “ 2 ( 2m ) ’ 


(7.23) 


(7.24) 


and therefore 

7 = 1 = 1-4, / = 1.90. 

At the other limit when the vibrational mode is fully excited, 


7.3?T. 7.3?. 9.3?. 

61 ~ 2 2m Cvl ~ 2 2m ^ “ 2 ( 2m ) ’ 


and hence 


Therefore, in general, 


7 = | = 1-29, / = 1.68. 

1.68 < / < 1.90. 


(7.25) 


In passing, we note that for monatomic gases, there are only three translational degrees of 
freedom and therefore 

_ 1 3?T 33?T 

62 2 m 2m ' 

Comparing (7.26) with (7.22) then shows that monatomic and diatomic gases have internal 
energies differing only on their zero point energies and that the basic assumption c v \ = c v 2 for 
an ideal dissociating gas is indeed valid. Also, the / value 2.5 for a monatomic gas can actually 
be read off from Eucken’s formula (7.21) by noting that c v 2 = 33?/(2m), c p 2 = 53?/(2m) and 
7 = 5/3. 
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We now return to (7.19). From the above discussion, it is appropriate to put c„ = 
33?/(2m), /i = 7/4, / 2 = 5/2. The coefficient of thermal conduction k for the mixture 

is then given by the formula 


k = 


213? 

8 m 


{i 


Pi 


+£. 


P2 


, '7 i + 2 * 1 . ja. [ • 

=1 M12 T *J 1 


(7.27) 


Specific heats of the mixture 


The specific internal energy for A and j 4 2 both have the same expression 

33?T 


e = 


2m 


The total internal energy of the mixture is then given by 

33?T\ 


E = 


2m 


(2m • rii + m • 712 ). 


(7.28) 


Since A and Aj do not necessarily have the same zero-point energy, relation (7.28) should be 
modified to 

3»T 1 

(7.29) 


_ 33?T. . 1 

E = ~2^-(2m • m + m • n 2 ) + -3?7dn 2 , 


where the last term represents the difference of zero-point energies of A over j4 2 . 

On dividing (7.29) by the total mass of the mixture (2mni +mn 2 ), we obtain the expression 
for the specific internal energy of the mixture: 


e = 33?T+ 3?r d a. 


(7.30) 


With the aid of this relation and the equation of state (7.6), we can easily calculate the 
specific enthalpy as follows: 


h = e + £ = (4 + q)KT + ftT d a. 


(7.31) 


To calculate specific heats, we have to first of all evaluate ( da/dT) p and ( da/dT) p . From 
(7-7), 

6>a. _ T d a(l - a) da _i.Mj.Hfh (7 

( dT' p T 2 2 -a ’ t' 1+ T' 2 ^ ^ ^ 7 ' 32 ^ 

It then follows from (7.30) and (7.31) that 


, de ^ , < o/ r <^2 «(!-«) 


% 

2' 


, Cp = (4 + Q )S + =(l + ^) 2 .a(l-a I ). (7.33) 


The quantity ( dh/dp)T which appears in (2.8) can now be calculated with the use of (7.31) 
and (7.2). The result is 


,dh s 5ft(T -f T d ) „ 2n 1 „ , „ 

( s? )r — •“(!-“ > = -2^(1 + y) •“(!-“)• 


(7.34) 
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7.2 Modification of the basic state 

As we have remarked before, the basic state equations (2.13) and (2.14) are independent of 
constitutive assumptions and therefore also valid for the gas mixture under consideration, 
although now p, p,, k and Cp are calculated from the more complicated expressions (7.6), 
(7.10), (7.27) and (7.33), respectively. In these expressions, the function a is given by (7.7). 
After non-dimensionalization, these equations become 



1 

P ~ (1 + a)T’ 

(7.35) 


(l + 77ll)T 3 / 2 i 2 (l + 77ll)f 3 / 2 

;i T + 77l! + -J- . 3 ^-(T + 7713 ) f + Th 2 + ^-~(f + m 3 )’ 

(7.36) 


T (l + mx)f 3 / 2 (10i 2 /7)(l + mi)f 3 / 2 

f + + + m 3 ) f + 7712 + ^— (? + rh 3 )’ 

(7.37) 


^ = 1 + ia + i(l + 3^)2 . a (i _ q 2 ), 

(7.38) 


a 2 pd ,= 

r = — Te 

l- a 2 poo 

(7.39) 

where 

S\ S 2 VSiSi 

m l = rp . m 2 = rp » m 3 = rp » 

loo loo loo 

- c * A-* - — - 
2 ~ Ai~ V2c 2 ’ 3 ~ Ai~V2c 2 2 

(7.40) 

For the purpose of asymptotic analysis, it is convenient to define two new constants a and b 

by 

Pd 2 l T d 

a ~ Poo M ' TooM 2 ’ 

(7.41) 

Then equation (7.39) becomes 

a 2 aT 

— p T 

1 — a 2 M* 

(7.42) 

which display 

s the physical fact that dissociation will take place in the hottest 

region where 


f = 0(M 2 ). 

An asymptotic analysis of equations (2.13) and (2.14) shows that the boundary layer struc- 
ture in the present case is similar to that for an ideal gas. In particular, the boundary layer 
can be divided into two regions: an inner region 77 = 0 (M -1 / 2 ) and a temperature adjustment 
region 77 = 0(1). In the inner region, we define new variables Y, T as in (3.5), but now (3.6) 
and (3.7) are replaced by 

( 1 +, 5 * l) (i 1 ^-S) ,+FF " = 0 ' (7 - 43) 
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where both D and fi are constants. Therefore, equations (7.43) and (7.44) are to be solved 
subject to the boundary conditions (3.8) and the asymptotic conditions (7.51). Results from 
such a numerical integration are shown in Fig.l together with those results corresponding to 
the undissociated model. 

In the temperature adjustment layer r\ = 0(1), a is exponentially small which means that 
no dissociation takes place in this region. The basic state equations (2.13) and (2.14) reduce 
to (3.13) and (3.14) which are appropriate for an ideal gas and which are to be solved subject 
to the same matching conditions (3.15) and the conditions (3.16) at infinity. 


7.3 Modification of the stability properties 

Since the boundary layer structure is similar to that for an ideal gas when dissociation is taken 
into account, the qualitative stability properties are also similar with appropriate quantitative 
modifications. In particular, we can still show that the mode trapped in the temperature ad- 
justment layer is most susceptible to Gortler vortices and the neutral Gortler number expands 
as in (3.35), but now the coefficient B is modified since it is the leading order contribution to 
the integration of basic state temperature across the whole boundary layer. As for higher order 
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correction terms to the Gortler number expansion, The effect of gas dissociation is dependent 
on the wall curvature. In the special curvature case, since the second term G in (4.1) is de- 
termined by solving partial differential equations in the undissociated temperature adjustment 
layer, it is not affected by consideration of gas dissociation; whilst in the more general curva- 
ture case, we can see from (5.9) that gas dissociation affects the second order correction term 
through E, but in the higher wavenumber case discussed in §5.3, the second order correction 
term is not affected by gas dissociation, as is clear from (5.44). 

In Fig. 19, we have plotted the variation of B with respect to the cooling coefficient both for 
the real gas model discussed here and for the ideal gas model used in section 3. It is clear from 
Fig.19 that the values of B is decreased by gas dissociation as well as by wall cooling. Therefore, 
both gas dissociation and wall-cooling are destabilizing. In our numerical integration of the 
boundary layer equations (7.43) and (7.44), we have taken rhi = rh= 0.508, A 2 = A 3 = 1, a = 
1.225 X 10 9 , b = 3.2653. We have also repeated our calculation for a few sets of different values 
for the above five constants. We find that the above prediction is still valid. 

In order to determine the effects of gas dissociation on the wall mode, we have integrated 
the perturbation equations (6.12)— (6.14) with the basic state given by the solutions of (7.49) 
and (7.50), subject to the same conditions as those used to produce the neutral curves in 
Fig. 13. Fig.20 to Fig.22 give a comparison of the neutral curves corresponding to the ideal gas 
model and the real gas model with dissociation. We observe that in each of these Figures, the 
two neutral curves corresponding to the two models intersect. Therefore, gas dissociation can 
have either a destabilizing effect or a destabilizing effect on the wall mode. 

8 Further discussion 

We have seen above that the Gortler mechanism in a hypersonic boundary layer of a Sutherland 
law fluid behaves in a predominantly viscous or inviscid manner depending on whether or not 
the wall curvature varies like where x denotes distance along the wall. When the wall 
curvature does not have this special form the vortices evolve over almost the whole of the 
wavenumber space in a nonparallel manner subject to viscous effects. The only exception to 
this case is at extremely high wavenumbers where the vortices evolve in a quasi-parallel manner 
essentially identical to that described for incompressible flows by Hall(1982). Since the special 
curvature distribution is possibly of little physical relevance we shall now concentrate on the 
results we have found for the more general curvature situation. 

The results we have found for the different wavenumber regimes for the general curvature 
distribution are summarized in Table 3 below. 
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Table 3 


wavenumber 

wall mode 0 ( M~ f ) 

strongly unstable inviscid mode 0(1) 
near neutral inviscid mode 0(1) 

nonparallel viscous mode 0(M? ) 

parallel viscous mode 0(M?) 


growth rate 
0 ( 1 ) 
0(MT) 
0(M?) 
0(M5) 
0(Mi) 


neutral value of G 

0(Mi) H 

no neutral G 
Gn + 0(M) + 

Gn + 0(M) H 

0(M2) H 


where Gjv is as defined below (3.35). 


The wall mode which was discussed in section 6 is in fact the counterpart of the so-called 
‘acoustic mode’ of the inviscid hypersonic instability theory of flat plate boundary layers. The 
latter Rayleigh instability has been discussed by Cowley and Hall(l990) and Smith and Brown 
(1990) for Chapman law fluids, and a limited discussion of this mode for Sutherland law fluids 
can be found in Blackaby, Cowley and Hall(1990). The acoustic mode and the wall mode 
discussed in section 6 have the property that they are concentrated in the wall layer where the 
streamwise velocity component of the basic state varies from zero at the wall to almost it’s free- 
stream value. The Rayleigh acoustic mode in general evolves in a quasi-parallel manner though 
in the presence of strong shocks this is not necessarily the case. The wall mode discussed in 
this paper evolves in a nonparallel manner and becomes progressively concentrated towards 
the edge of the wall layer as the local wavenumber increases. 

When the wavenumber of the vortex becomes 0(1) then the disturbance modifies itself so 
as to become concentrated in the adjustment layer where the basic state temperature adjusts 
rapidly to it’s free-stream value. The counterpart of this mode in the Rayleigh instability 
problem is the so-called ’vorticity’ mode investigated by Blackaby, Cowley and Hall (1990) for 
Sutherland law fluids, and by Smith and Brown (1990) for Chapman law fluids. We found in 
the present paper that when the Gortler number G is as given by (3.27) with G* > Q and Q 
given by (3.28), then the appropriate expansion of the disturbance field is given by (3.31). Thus 
the mode has growth rate 0(M*) and the growth rate is given by the solution of the eigenvalue 
problem specified by (3.32) subject to the condition that Vo should vanish at 0,oo. Fig. 3 
shows that the growth rate increases monotonically from zero as the wavenumber increases 
and tends to a constant at large wavenumbers. At large wavenumbers the growth rate can be 
matched onto the parallel viscous mode growth rate when the vortex wavenumber becomes 
formally 0(M t). We note here that the inviscid mode matches directly onto the parallel 
viscous mode structure at sufficiently high Gortler numbers. This is exactly the situation with 
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the inviscid- viscous connection between temporally growing Gortler vortices in incompressible 
flows and is a direct consequence of the fact that the growth rate of the inviscid mode shown 
in Fig.3 tends to a constant at large values of the wavenumber. In the corresponding spatially 
growing problem for incompressible flows Denier, Hall and Seddougui (1990) show that an 
intermediate region is required to make the inviscid-viscous connection and indeed that the 
maximum possible growth rate is achieved in that interval. We stress that this intermediate 
regime has no connection with the nonparallel viscous mode discussed in section 5. The latter 
regime is appropriate only to the inviscid-viscous connection problem at mildly supercritical 
values of G. In Fig. 23 we have sketched the dependence of the growth rate on the the 
wavenumber for the cases when the Gortler number is mildly supercritical (ie when G differs 
only slightly from its value required to overcome the strong 0(M*) curvature of the basic 
state) and the strongly unstable case with G ~ 0(M a). The broken parts of the curve denote 
regimes where the growth rate is given by a nonparallel calculation. In Fig.24 we sketch the 
neutral curve in the local Gortler number-local wavenumber plane. In this figure we have 
indicated the corrections to the neutral curve associated with nonparallel effects. 

There are no available experimental results with which we can compare with our calcula- 
tions, this is because experiments at the high Mach numbers, say 10 — 30, appropriate to our 
work are exceedingly difficult to perform. Thus for design purposes it is fair to say that a 
theoretical approach is the only means at this stage to predict the likely evolution of Gortler 
vortices in growing hypersonic boundary layers. Current transition prediction methods are all 
based on some amplitude growth criterion based on the linear growth of a disturbance. We 
have seen above that for a realistic hypersonic boundary layer the regime where vortex growth 
is likely to occur is the strongly unstable inviscid one. Hence in any transition prediction 
method for a hypersonic boundary layer it would be appropriate to simply compute the local 
Gortler number and obtain the corresponding growth rate from Fig.3; this of course should 
only be done for wavenumbers less than the neutral value associated with parallel viscous mode 
structure. 

Finally we close by making a few remarks about the results of section 7 which concerned the 
effect of gas dissociation and wall cooling on vortex growth. First let us make a few remarks 
about the effect of wall cooling, since the most unstable vortices correspond to the strongly 
unstable inviscid mode we concentrate on that situation. In that case the main effect of wall 
cooling or real gas effects is to alter the quantity B which fixes the scaled Gortler number 
above which instability can occur. Fig. 19 shows that for an ideal gas the effect of wall cooling 
is to reduce B by a factor of about .5 when the wall temperature is cooled from it’s adiabatic 
value by a factor of 10. We find that the eigenvalue problem associated with (3.32) has the 
maximum value of (3 altered only by a small amount when this happens so that at a given 
value of the Gortler number the growth rate is increased by a factor of about 1.4 when the 
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wall temperature is reduced by a factor of 10 from it’s adiabatic value. By contrast real gas 
effects have a negligible impact on the growth rate of the strongly unstable mode. Thus for 
example in Fig. 19 we see that, at a fixed value of the cooling coefficient, B varies by only 
about 10 percent when the real gas model is used. Thus the critical Gortler number or the 
disturbance growth rate is altered only slightly by real gas effects so it would seem that any 
transition prediction method could quite sensibly ignore such a complication; certainly the 
error associated with doing so would be negligible compared with the inherent error of the 
prediction method. 
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Appendix: Basic state properties 


Listed here are various derivatives related to the basic state, which are used frequently in the 
determination of the relative orders of the perturbation quantities. 

The temperature adjustment layer 
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where 


u* = F"(Y)Y X , Uy = F"{Y)Y V , 

(A9) 

M 3 / 2 - T'F F" „ 

% = - { 2x)z/2 [ - TF + iF '- f )0+ r n] ’ 

(A10) 

v y = ^[-f' F +F"D], 

(All) 

_d _d , d f d 

U dx +V dy ~ F dX 2 xdY' 

(A12) 

M 3 / 2 

uv x + vv y = -j^-^[-TFF' - T'F 2 + (F 12 + FF")Q], 

(A13) 

UU X + VUy = , 

(A14) 

Q(Y) d = [ Y TdY. 
Jo 

(A15) 


REFERENCES 

1. Becker E. (1968) Gas Dynamics. Academic Press, New York, London. 


2. Blackaby, N. D. ,Hall, P. and Cowley, S. J. (1990). On the instability of hypersonic flow 
past a flat plate. ICASE Report 90-40 and submitted to J. Fluid Mech. 


3. Buddenberg J.W. and Wilke J.R. (1949). Calculation of gas mixture viscosities. Indus- 
trial and Engineering Chemistry 41, 1345-1347. 


4. Chapman S. and Cowling ,T.G. (1970), The Mathematical Theory of Non-Uniform Gases, 
3rd ed. Cambridge Univ. Press, London and New York. 

5. Cowley, S. and Hall, P. (1990). On the instability of hypersonic flow past a wedge. J. 
Fluid Mech. 214, 17. 

6. Denier, J.P., Hall, P. and Seddougui, S. (1990) , On the receptivity problem for Gortler 
vortices: vortex motions induced by wall roughness. ICASE Report No. 90-31. and Phil. 
Trans. Roy. Soc. A, In Press. 

7. Freeman, N. C. and Lam, S. H. (1959) , On the Mach number independence principle for 
a hypersonic boundary layer. Princeton University Report 471. 

8. Hall P. (1982). Taylor- Gortler vortices in fully developed or boundary layer flows: linear 
theory. J. Fluid Mech., 124, 475-494. 

9. Hall P. (1983). The linear development of Gortler vortices in growing boundary layers. 
J. Fluid Mech. 130, 41-57. 


47 



10. Hall, P. (1988). The nonlinear development of Gortler vortices in growing boundary 
layers. 

11. Hall, P. (1990) Gortler vortices in growing boundary layers :the leading edge receptivity 
problem, linear growth and the nonlinear breakdown stage. Mathematika, In press, (based 
on ICASE Report 89-81) 

12. Hall, P. and Malik, M.R. (1987). Gortler vortices in compressible boundary layers. J. 
Engrg. Math. 23, 239 

13. Lighthill M.J. (1957), Dynamics of a dissociating gas: Part I Equilibrium flow. J. Fluid 
Mech., 2 , 1-32. 

14. Lindsay A.L. and Bromley, L.A. (1950). Thermal conductivity of gas mixtures. Industrial 
and Engineering Chemistry 42, 1508-1511. 

15. Smith, F. T. and Brown, S. N. (1990). The inviscid instability of a Blasius boundary 
layer at large values of the Mach number. J. Fluid Mech. 219, 499. 

16. Spall, R. E. and Malik, M. R. (1989). Goertler vortices in supersonic and hypersonic 
boundary layers. Phys. Fluids A 1, 1822. 

17. Stewartson K. (1964). The Theory of Laminar Boundary Layers in Compressible Flows. 
Clarendon Press, Oxford. 

18. Wadey, P. (1990). Gortler vortices in compressible boundary layers. Exeter Univ. PhD 
thesis and submitted for publication. 


48 




Fig.1 Basic state temperature in the wall layer corresponding to 
two models. Here n is the cooling coefficient. 
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Fig. 2. Basic state temperature in the adjustment layer. 
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Fig. 3 The dependence of the eigenvalue y/2x(3 2 /(G* - Q ) of (3.32) on the 
local wavenumber fc, which agrees with the results from asymptotic analysis 
that as k — > 0, \/2 xf3 2 /(G* -Q) is of order k; whilst as k -* oo, y/2 x(3 2 /(G" - 
Q) tends to the constant 0.5786. 
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Fig. 4 The eigenfunctions of (3.32) corresponding to a set of decreasing 
wavenumbers k. Asymptotic analysis shows that as k -* 0, vortices are 
confined to a wall layer of 0(k x l 3 ) thickness and Vo tends to a constant 
when moving away from the wall layer. 
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Fig. 5 The eigenfunctions of (3.32) corresponding to a set of increasing 
wavenumbers k. Asymptotic analysis shows that as k — ► oo, vortices are 
confined to an 0(k ~ */ 2 ) thin layer centred at i) = 2.3228 where —T'/T 2 
attains its maximum. 


53 






Local Gortier number 


— • E5 as energy measure 
• • E4 as energy measure 
— E3 as energy measure 
— — E2 as energy measure 
•••• El as energy measure 
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Fig.6. Neutral curves obtained when we impose the following initial coni- 
tions at the upstream postion Xo = 60: 

U = e- 1 / T)4 • e-^/r?, V = 0, 6=r)U. 





Local Gortler number 


4000 



Fig. 7. Neutral curves obtained when we impose the following three initial 
conditions at the downstream location .Yo = 20 and use E$ as the energy 
measure to monitor the energy growth. 

(i) . U = e" 1 /” 4 • V = 0, 6 = tjU. 

(ii) . U = ^e" 10 ^ 5 ■ e-”, V = e~ w '^ ■ e~\ 9 = 0, 

(iii) . U = e~ 40 ^ 6 • e-», V = i fe- 40 '"* ■ e~'\ 0 = 0. 
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Fig. 8 . Neutral curves obtained when we impose the following initial con- 
ditions at three different upstream positions A'o = 20, 40, 60 and use E$ as 


energy measure: 

U = e- 1 ^' -e-'/r), V = 0, § = t)U. 
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Fig. 12 Temperature perturbation at different downstream locations 
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Fig. 13 The neutral curves for the wall mode. 
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Fig.20 A comparison of neutral curves corresponding to two models 

(El as energy measure) 
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Fig. 21 A comparison of neutral curves corresponding to two models 

(E2 as energy measure) 
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Local wavenumber 

Fig.22 A comparison of neutral curves corresponding to two models 
(E3 as energy measure) 
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Fig. 24 The neutral curve for Gortler vortices in hypersonic boundary layers. The dashed part of the curve 
indicates that at the given wavenumber the curve depends at zeroth order on nonparalle effects. The verical 
arrows indicate the size. of the next order correction to the neutral curve, apart torn the regime where the 
wavenumber is 0[M‘) the correction depends on nonparallel effects. 
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